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SYNOPSIS 


™ MASS-POIITT IICSBL 

Gupta ( 1969 ) obtained catisdactory values of E(I)/®(2), 

Ij = 4j6, etc., foy tlio deformed nuclei in 'zhe rrosion 150 < A < 1! 
by representing a typical nucleus by a twc mass-point system 
with, a harmonic restoring force a.s a function of the distance 
between them and using semi-classical treatment. Although the 
model is not applicable to spherical nuclei, for deformed nuclei 
in the rare-earth region, it correlates the que/.rupole moment 
and the moment of Inertia, data more satisfa,ctorily than the rigii 
body or the hydrodynamic models. Such a representation of a 
nueleus is expected to be more appropria.te for very large 
deformations, when the nucleus is about to undergo binary fissioa 
However, what is rather interesting to note is tli:,t it is also 
applicable to nuclei v/ith relatively sma.ll deformations. 

Trainor and Gupta (1971) modified the above- model to a 
governor neidel, so as to obtain zero moment of inertia for 
spherical nuclei, assuming that e. central spherical portion of 
radius equal to the semi— minor eacis of the spheroid c^oes not 
participate In the rotational motion. The model was found to : 
give satisfactory results. 

In view of the above, it was felt that a quantum ^ 

mechanioal treatment of rotational-vihrational motion of the 
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two mass-point system with harmonic and with anhamonio 
potentials woxild yield better results 0,0 composed to the results 
of Gupta (1969). The energy expression thus obto.ined v/as found 
to contain in addition to the seni-classica,l terms (Gupta 1969), 
a beta vibrational energy term sjid 0 , 1-dependent zero point 
vibrational energy term, both contributing significcaitly to the 
energy of high angulo,r momentuim roto.tional states. Several 
quantitif(e calculated on the basis of the model v;ere found to be 
in better agreement with the corresponding ergperimental values. 
Besides, st. number of features characteristic to these nuclei as 
obtained by Hariscotti et al (1969) and Mosel and Greiner (1968) 
are also observed for the tv/o mass-point model. The vo.riation 
o^^ the stiffness with the proton ojid the neutron numbers is of 
particular interest (Prakash et 0,1 1973). 

Further, we hate re-examined the governor model of Traino] 
and Gupta (^971), by choosing tl\e radius of the rotationally 
invariant core as a variableu such that the two mass-point 
model and the governor model are obtained with appropriaite 

f 

limiting values of the radius of the core. The varir.tion of 
relative densities inside and outside the spherical core as a 
function of its radius is rather interesting, 

! 

2 . ROTATIQIAI FLOW MOIBL ' 

As an improvement over the two mass-point model, we 
reconsidered the continuum model as proposed initially by Bohr 


XTl 


and Mottelson (1955). Tlie flovr of the continu-uon is laiovm to he 
neither irrotational nor rigid. 7fe, therefore, coiisidered a 
rotational flow model (RFK) with a finite hut unifonn vorticity 
throughout the nuclear -volume. The treatment wo.s carried out 
in two dimensions with a subsequent generalization to flows in 
three dimensions. The velocity field for the case of pure 
rota.tion of the ellipsoidal boundary wa-s uniquely found. 


With the knov/ledge of the flow, the Hamiltonian of the 


system could he written down. It was sho\m that the Hamiltonian 


reduces to a rotational form H = E (lf/2Ip^) only if the 


D=1 


Ej 


vorticity components c. are proportional to the rotational 

J 


frequency w. of the ellipsoido,! boundary. It is thus concluded 
J 


that in the continuum model, the intrinsic motion is strongly 
coupled v.dth the rotational motion, and the adiaha,tic approxi- 


mation is valid only if the proportionality constants (Ic . ) are 

J 


much larger than unity. Ass-uming c- = Ic- w. an expression for 

d J J 


the effective moment of inertia I-^^ as a function of k. and 

d 


the shape parameter = h^ ms obtained, Eor k^ = 0 and 1 , 

the flow corresponds to the rigid and the irrotational flow 
respectively, and it is found that the effective moment of 
inertia also goes over to the appropriate limiting values. 


It was found tha.t the experimental values of thb moments 
of inertia and ths quadrupole moments in the ground state 
rotational band of deformed even-even nuclei in the rare-earth 



region could be accounted for by clioosing the vali’.eo the 

vorticity constant (k-) in the rcjige 0.8 to 1.0. It was found 

J 

that the va^lue of the vorticitj?' constant goes to 1 , giving 
irrotational flow, as the nuclear deforma.tion goes to zero. 

Thus, it is concluded that 'che nuclear flow v/liic'i is rotational 
for deformed nuclei tends to become irrotatioiial for spherical 
nuclei, with the consequent result that spherico.l nuclei do not 
show rotational spectrum. 

The hi^er rotational bands have been analysed in terms 
of processional motion of the ellipsoidal boimdar 3 ’' \;ith rotationa] 
flow. Small but finite values of moment of inortic, a,bout the 
longest axis in the ellipsoidal drop has been accotmted for 
'(Prakash et al. 1972, Prakash anri Pe,-h panda 197'.). 

Attempt is made to obtain a univercn.1 re3-ationship betwe^, 
the vorticity and the shape paremoter b. valid for a.ll rota- 

J 

tional bands. Attempt is also made to write the liamiltonian as 
a function of deformation from the determined flow pattern in 
order to estimate enuilibrium deformation and vibrational stiffnes 
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CHAPTER-I 


INTROIiUaTIOR 

Kienonieiiological models liave played a useful role in 
understanding nuclear structure ever since the introduction 
of the liquid drop model hy Hiels Bohr and Kalcker (1937). 

With the subsequent development of the single particle model 
independently by Mayer (1949) and Haxel» Jensen ana Suess 
(1949), foundation was laid for the microscopic trea,tment of 
the nuclear system. While this model was eminently successful 
in accounting for the observed shell effects > processes like 
the fission continued to be studied in terms of the continuum 
model. While the understanding of the independent particle 
motion in terms of intemucleon forces improved with the advent 
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of Bruecloiei*’ s flieory (1954) » i't was also 'becoming apparent 
that tbe nuclear system was capable of collective motion 
(Bohr 1952, Bohr and Mottelson 1953)» such as rotations and 
vibrations. The relatively simple and common features of 
collective motion of a large ntmiber of nuclei permitted appli- 
cation of phenomenological models with advantage (Sood 1968). 

In such treatments it was not intended tliat the complicated 
motion of the nucleons be described, Q3aey were intended for 
classification and for inter-relating the large wealth of 
collective data. Such analysis brings out broad trends which 
are of interest in themselves and which provide guidelines for 
the microscopic treatments. In the recent past, significant 
advancement has been made in accountirg for the collective 
behaviour of the nuclear system in tei ms of microscopic theories, 
while phenomenological treatments continued (Diamond et al 1964f 
Krutov 1968, Mariscotti et al 1969, &upta 1969 and Irainor and 
G-upta 1971) to bring out inter-relating features, 

!Ehe present v/ork consists of exploration of some pheno- 
menological macroscopic models to account for the rotational 
spectra of even-even nuclei primarily in the rare-earth region. 
In this region with 150 < A < 190, rotational levels with 
spacing of the form L(L+1) have been observed (Lieder et al 
1971) up to I * 22. Ihe simplest model which could be consi- 
dered here is that of a rigid rotator. The broad features of 



the ohserred spectra, such as L(L+1 ) dependence, and the 
existence of the higher precessional hands can indeed he 
accounted for, on the basis of this model. However, there is 
evidence that nuclei are not rigid, fhey can vibrate, and 
their moment of inertia increases (Sood 1968) v/itli the spin I, 
so that there is noticeable departure from 1(1+1) rule, as 
in the case of molecular rotation. The experimental value of 
the moment of inertia as obtained from absolute level spacing 
is found to be significantly smaller than the rigid body value 
obtained from nuclear mass and size, finally, it is Icnown 
that only the nuclei with finite deformation exhibit rotational 
spectra which is not true of rigid rotators. 

Bohr and Mottelson (1953) considered the nucleus as a 
continuum, characterised by flows. Rigid rotation is then 
one possible flow of such a continuum. They asstimed that, the 
continuum is incompressible and that the flow is irrotational. 
The assumption of the irrotationality leads to possible shape 
vibrations but no rotation for nuclei with spherical equilibrium 
shape, in accordance with experimental observations. The Bohr- 
Mottelson treatment, however, led to vibrational frequencies 
which are too high and to values of moment of inertia which are 
too low. The experimental values of the moment of inertia, 
obtained from the observed rotational levels, therefore, lie 
between the rigid body values and the irrotational. flow values. 
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Miile the Bolir-Iiottelson continum model over-estimated 
the absolute spacing, it provided a basis for the consideration 
of finer details of the relative spacing, IDhus, the treatmeni; 
assuming spheroidal shapes has been extended (Pal 1971) to 
ellipsoidal shapes. The increa,se in the moment of inertia with 
the angular momentum has been considered (Sood 196?, Gupta 1967, 
1969) in terms of the centrifugal stretching of the non-rigid 
continuum. 

Other phenomenological treatments have also been made. 
Thus, satisfactory fits have been obtained (Mariscotti et al 
1969) by treating the moment of inertia as a variable parameter. 
Gupta (1968) has considered the deformed nucleus in terms of 
a two mass-point system. Trainor and Gupta (1971) have 
considered a governor model, vhere a d-^formed nucleus is con- 
sidered as consisting of a non-rotational spherical core and 
the extra- spherical region responsible for the rotational 
moment of inertia. 

The present work consists of two parts. In the first 
part, we have applied the two mass-point model to the even-even 
nuclei in the rare-earth region and have re-examined the governor 
model. In the second, we have considered the continuum model 
with flows of finite vorticity. The vorticity is assumed to be 
finite but uniform tliroughout the nuclear volume. The flow for 
the case of pure rotations is determined. The coupling of 
intrinsic flow (in the body fixed frame) with the rotational 



motion, is examined, An expression for the effective moment of 
inertia is obtained. Comparison with experimental results shows 
that the flow is rotational for deformed nuclei aiid it goes over 
to irrotational flow for spherical nuclei, Highel* rotational 
hands have been analysed to account for the small but finite 
values of the moment of inertia about the longest axis of the 
ellipsoid. Attempt is made to obtain a relationship between 
the vorticity and the deformation and thus to write the Hamilto- 
nian in terms of the deformation in order to estimate the equi- 
librium deformation and the vibrational stiffness. 

In Chapter II, the available experimental data and the 
relevant theoretical methods to study the collective motion in 
even-even nuclei in the rare-earth region is discussed, • 

In Chapter III, ihe main features of the rigid body motion 
with its success and failures as applied to nuclear collective 
motion is presented. 

In Chapter 17, the two mass-point representation of a 
defoamied even-even nucleus is considered quantum mechanically. 

In Chapter 7, the general features of the continuum model 
are outlined. 

Chapter 71 presents the hydro dynamical model with a 
rotational flow and the last Chapter contains a summary of 
results* 
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CHAPTER - II 

OOLLEGTIYE MOTIOH; E^ERIMEITTAL ARE !IHBORETIQAL STATUS 

Tlie primary criterion for iiie success of a model proposed 
to account for the collective phenomenon in even— even nuclei in 
the rare-earth region, is to account satisfactorily for the 
available experimental data and their systematics. Hence, it 
will be appropriate, here, to discuss briefly some aspects of 
modern experimental methods employed in such studies, the 
systematics of the available experimental data and the theoretical 
methods to explain them. 

11,1 Experimental Methods 

The most promising tool^ till date, for the exploration 
of the nuclear rotational and vibrational states is the Coulomb 
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excitation of a target nucleus iDy the electromagnetic field of 
charged projectiles which oxe not capable of penetrating into 
the region of nuclear forces. The information obtained out of 
such experiments regarding various nuclear properties such as 
spins and parities, transition rates, guadrupole moments, 
moments of inertia, magnetic moments, etc., are very accurate, 
for the theory of Coulomb excitation is accurately understood 
within the framework of electro d 3 ?namics* A collection of 
important papers and review articles on this subject is given 
by Alder and Winther (1966). 

Recently, the use of heavy ion as a projectile has proved 
to be most usefixl in the study of collective phenomenon in nuclei. 
In this case, the electric field exerted by the projectile on 
the target is so large that higher order excitations occur. As 
the heavy ions carry large orbital angular momentum of the order 
of 50h to 100h, excitation of rotational states of angular 
momentum as high as 20h has been possible. The process often 
results in multiple Coulomb excitations of the target nucleus, 
wherein the successive levels are excited by its predecessor, 
or eventually result in very hi^ excitation of the target nucleus 
and then proceed to the evaporation of several neutrons or charged 
particles. This method has been very useful in the production 
and study of neutron deficient nuclei and has contributed immensely 
to the wealth of experimental data of such nuclei in the region 
150 < A < 190 , An excellent review of heavy ion reaction has been 
given by Newton (1970). 
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]?urther, it has heen foTmd that the cross-sections and 
ganma ray angular distributions in Ooulomb excitations depend 
to a measurable degree on the static quadrupole moments of 
nuclear levels. This is called the ’re-orientation effect’ and 
has recently become an important tool for measuring the static 
quadrupole moment^ of excited levels, especially, of some ’one- 
phonon’ states in vibrational nuclei. A good review of this 
subject has been given by de Boer and Bichler (1968)* 

II. 2 SHJBY OB E3CPERIMMlAIr BAIA 

(Che experimental methods outlined above and many others 
finally provide us with the following experimental infoimation: 

(a) spins and parities of various levels 

(b) electro-magnetic transition rates 

(c) static and d3mamic quadrupole moments 

(d) magnetic moments 

(e) energies of various states and corresponding moments of inertia 

11,2,1 Spins and Parities 

In all even-even nuclei, following Pauli exclusion 
principle, the even number of protons and even number of neutrons 
will separately pair off among themselves to give a zero contri- 
bution to the total spin of ihe nucleus with a positive parity 
in its ground state. The ground state with l’*' = o'*’ forms the 
head of a set of levels with seq.uence 1^ = 2 , 4 , 6 , etc,, vhich 
constitute the ground state rotational band. Each time a spin 
repeats, it forms the head of a new excited band with a set of 
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levels luilt upon it with = 2"^. At sufficiently hi^ 

excitations ('v 1 MeY) a set of levels with 1^ = 1*^, 2'*’, 3^3 
etc,, have also been observed, 1^ = 1"^ forms the head of yet 
another excited band and so on. The spin of the level forming 
the head of a particular band is usually denoted by K to 
distinguish it from the rest of its members. It will be shown 
later that K is nothing but the projection of Ii along the 
symmetry axis of a nucleus and is always less than or equal to I, 

11.2,2 Electro-magnetic transition Rates 

The nuclei in the rare-earth region exhibit electric 
quadrupole transition rates which are about 100 times those of 
single particle estimates and hence is attributed to the 
collective excitation of a large number of micleons. The 
reduced S2 transition probabilities j’or rotational excitations 
in even-even nuclei of spheroidal shape are given by (Bohr and 
Mottelson 1955) 

B(E2) = -If - aS ■(I-+1-S.) (l+H-K) iL+2-K) U.+.Z+K) 

° (Ii+1) (21+3) (1+2) (21+5) 

for transitions of the type L+2 ->■ B for a given K, Here 
Qq is the intrinsic (static) quadrupole moment in the state L. 
Eor L = K = 0, eq.(l) gives a relationship between the B(E2) 
and Qq for the ground state of an even-even nucleus. In 
fact, B(E2) is related to the partial E2 gamma ray half-life 
(T^y- 2 ) and the transition energy Ey by the expression (L*<5bner 
et. al. 1970 ) 
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B(E2) = 56.56/(E^T^/2) 

for transitions of the type 2*^ -* o'*". Here B(E2) is given 
in e^ 10“^® cm^, seconds, and Ey in KeV, 

Eq_. (1) can he used to compare B(E2) values between 
any two pairs of levels to obtain values of branching ratios. 

11*2*3 Quadrupole Moments 

The quadrupole moment of a nucleus arises out of the 
non-spherical distribution of cho,rge and hence is a measure 
of the degree of departure of the nucleus from sphericity. 

The quadrupole moment defined with respect to an axis fixed 
in the body of the nucleus is called the intrinsic or the 
static quadrupole moment (Q^) and that defined with respect to 
an axis fixed in space is called the dynamic or the spectroscopic 
quadrupole moment (q)* Hence, the value of Q for a given 
state is determined by the spin and the distribution of charge 
of the nucleus in that state, ?diereas Qq is determined only 
by the latter. In fact, for spheroidal nuclei, the two are 
related by the expression 

Thus, in the ground state of a deformed even-even nucleus 
(1 = K = 0) although Q is zero, is not. 

Assuming the nucleus to be a uniformly charged body of 
ellipsoidal shape v/ith a.| , a, 2 > O’® principal axes, the 
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following expression for is obtained (Kumar 1972) 


Qq/® 


(Z/5) ( 


2 2 

2a^ “ a^ “ n, 


|) + f6 (ag - a^) J 


Hie two terms vfithin the square hrackets respectively arise 
out of the axially symmetric deformation (p) and non~axially 
symmetric deformation (H ) of the nuclear surface 

i (x?/a?) = 1 

3=1 ^ ^ 

with a^ > ag > . Introducing polar coordinates and solving 

for R(9,(|)),. the parameters a^ , ellipsoid can 

he obtained in terms of the parameters p and (Kumar 1972). 

a^(p,Y) = Rq (p,'/ ) [ l+V^Ti p Cos (/ + 2rtD/3) (5) 

/6 

where Rq(P,'/ )= Rq ^ “ (15A7t)p^ + I* Cos^f] 


8,^ 8,r) Q**; 


Here R.^ is the radius of the equal volume sphere. The ellipsoidal 
surface defined hy eq. (5) is not identical to the quadrupole 
surface defined hy Bohr and Mottelson (1953). Ho?/ever, for 
small spheroidal deformation (V = O), the parameter p is the 
same as that defined hy Bohr and Mottelson, Assuming the 
nucleus to he a rigid prolate spheroid (a^ = ag and */ = 0), 
for small symmetric defoimation, eq, (3) combined v/ith the 
eq, (5) approximately reduces to 


(3/f^) ZR^ P(1+ 0.16 p) 


( 6 ) 
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Tile large valties of B(S2) for all even-even nuclei 
in tlie region 150 < A < 190 show that they are strongly 
deformed hoth in their ground and excited states. The values 
of p obtained from the loiown values of for the ground 

state range between 0.16 and 0.35} and is maximum for the 
nuclei in the middle of the region. The experimental data 
is available only for the ground state (0 ) and the first 
excited rotational state (2 ). The data for others is scarce. 

It ViTill be shown at the end of this section that there 
exists a definite correlation between the intrinsic quadrupole 
moment (Qq) and the energy of the corresponding state. An excellent 
compilation of recent data on Q^, obtained from various experiments 
has been given by IWbner et. al. (1970). 

H»2.4 Magnetic Moments 

The magnetic moment ) expressed in units of the 
nuclear magneton is 

= • 1 > / (L + 1) . (7) 

^ere ¥-° = 5 + I + gj^S . 

Here, g^ and g^ are respectively the g-f actors for the 
intrinsic spin s and orbital angular mcmentum 1 of the 
nucleon, and is that for the rotational angular momentum 

R of the nuclear surface. According to the above expression-, 
all the even-even nuclei must have ^ = 0 in their ground state 
since their ground state spin is zero. This is true for all 
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even-eYen nuclei in the region 150 < A < 190. 

3!he experimental data availahle for nuclei of our' int-erest 
is very sparse. In most of the determinations, the sign of 
is uncertain and only for three nuclei ^^^hy; "^^^Hy) the 

value of ^ is knom hoth in magnitude and sign. Hence, at 
present, it is hard to draw any conclusion regarding the 
systematics of yti, in excellent compilation' of the experimental 
data on obtained from various experimental methods has been 
given by Itiller and Cohen (1969). 

11,2.5 Systematics of Level Energies j 

We consider all the even-even nuclei lying between the 

closed shells H = 82, Z = 50 and H = 126, Z = 82. Por nuclei 

near the closed shells, the first excited state is a singlet 
« "I" 

with spin-parity 2 , the second is a triplet with spin parity 
0 , 2 , 4 J and it is experimentally found that the ratio of 
the energies between the two is about 2,2 thereby implying 
that these states may be vibrational in character. However, the 
degeneracy of the triplet state has been observed experimentally. 
As we move away from the closed shells, the energies of first 
and second excited states go progressively lower and the triplet 

N 

state splitting becomes more and more prominent. 

Nuclei in the region 158 < A < 186 are found to hhve a 

^ 

set of low-lying excited states with spin-parity 0 , 2 , 4 » 6 , 
etc,, with the energy of the lowest excited state (2 ) to be 
of the order of 100 KeV, The energies of these states, to a first 
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approximation, are foimd to "be given "by tlie familiar 
rule such that E(4'^)/E(2'^) = 10/3, B( 6‘^)/E(2'^) = 7, 

E(8'^)/B( 2"^) = 12_ and so on; implying that these states arise 
out of the rotations of the nucleus as a whole. However, for 
states with highex- L the experimental energies go progressively 
helow those expected from L(l+1) rule, indicating thereby 
that there is an increase in the moment of inertia of the 
nucleus, because a rotating nucleus can undergo centrifugal 
stretching. Usually, for a deformed nucleus the gromd state 
moment o-p ine'^'^'a is taken to be approximately equal to 

that calculated using the energy of 2”** state and the 1(1+1 ) 
rule, Eor lar^/'^ valxes of 1 the value of i^^pt ’ calculated, 

is foimd to approach the rigid body value 1 ^, 

A plot of the experimental values of quadrupole moments 

and energies of the 2 state as a function of neutron number 

(A-Z) for all even-even nuclei in the rare-earth region shown 

in EigureH.I indicates that there is a definite correlation 

between the two data. Also, we see that for a given set of 

isotopes (same Z) in the region 84<( A-Z)<102, the energy of 
+ 

2 state is found to increase vath the decrease in the neutron 
number v/hereas the corresponding quadrupole moment increases 
with the increase in the neutron number H, Eor nuclei in the 
region 102<(A-Z)<1 1 6 the reverse of this case is observed. 

Also, it is clear from Eigurelll that there is an asymmetry 
in some of the properties of the nuclei belonging to the two 
families. 
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Another rather interesting and important systematics 
recently foimd in these nuclei first by Johnson et, al. (1972) 
and subsequently by a number of others (Thieberger et. al. 1972; 
Lieder et, al, 1972) is that the -variation of moment of inertia 
(Ig;^^) as a function of the square of rotational frequency (-flw) 
for all the states in the ground state rotational band follo-ws 
a typical S-type curve for a given nucleus for which the energies 
of sufficiently high angular momentum states (li^ > H"*") are 
known, The two quantities plotted in figure 11,2 for some of the 
rare-earth nuclei are given by (Appendix I) 


2 2 
(Hw)^ = [e(L) - E(L~2)1 KeV‘ 


( 8 ) 


vhterc I can take the values 2,4,6, etc. 

In addition to the low-lying levels described above, 
often groups of excited states (at energies around 1 MeT) with 
spin-parity o"^ , 2"^ , 4"^ , etc., and o"^ , l"*" , 2 ^ , etc., have 
also been experimentally observed. However, the experimental 
data in this part of the spectrum is so scarce that it is hard 
to draw any conclusion regarding their systematics, A very 
good compilation of the recent experimental energy level data 
for various rotational bands is given by Mitsuo Sakai (1970), 

Thus, the spectrum exhibited by most nuclei near the 
closed shells is of vibrational type and as we go toward the 
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middle of the region there is a gradual transition from the 
vibrational type to the rotational type. However, the nuclei 
in the transition region, though a few in number, are found to 
exh_ibit both vibrational and rotational characteristics and 
hence are relatively complex to analyse. 

II. 3 IHEORETIOJLL STUDIES 

« Basically there are two different approaches generally 
adopted in the study of nuclear collective motion: (i) microscopic 
and (ii) macroscopic. 

The current microscopic investigations are mostly based 
on the shell model picture of the nucleus wherein it is assumed 
that each nucleon moves in a deformed effective field (this 
makes the model .someTiiiiat collective in nature) due to the rest, 
since the exact form of nucleon-nucleon interaction is not known, 
Althou^ the model very successfully explained the experimental 
data and their trends in case of light nuclei (A < 40), it failed 
to do so for the nuclei in the region 150 < A < 190 and A > 220. 

An excellent discussion of various microscopic theories of 
nuclear collective motion is given in the book by Rowe (1970), 

Macroscopic studies of nuclear collective motion are 
based on the liquid drop model of the nucleus proposed by 
Riels Bohr and Zalckar (1937) axid its subsequent developments. 
Perhaps, the most interesting aspect of its development has 
been Ats gradual merger with its chief competitor, the independent 



particle model giving rise to the unified model (Bohr 1952, 
Bohr and Mottelson 1953)? which could descrihe the collective 
and particle phenomena simultaneously. This idea has been 
developed considerably in recent years. In fact, the concept 
of particle interaction via the fields that they generate, 
which need not be static as in siii 5 )le shell model, underlies 
virtually the whole of modem collective theory. Thus, the 
macroscopic models, although phenomenological in nature, have 
been instrumental in the development of microscopic theories. 

As the work presented in this thesis is based on the 
phenomenological models of nuclear collective motion, it will 
be of use, here, to discuss the salient features of some of the 
important and recent phenomenological models. 


II. 4 STUDY OB MCROSCOPIG MODELS 

Por axially symmetric even~even nuclei in the rare-earth 
and actinide regions, the energies of rotational levels in the 
ground state band (K = O) are approximately given by 

E(L) = (*^21) D(D+1) (10) 


v&LQTQ L = 0,2,4, etc,, and I is rotational moment of inertia 
of the nucleus, the experimental value (l^^^) of which is 
obtained using the energy B(2) of the first excited 1=2 state. 


In quantum mechanics, the rotational moment of inertia of 
spherical systems is zero. In fact, since all directions in 
such a body are indistinguishable, it follows that it cannot be 



"brought into rotation (Uemirovskii I960), 

11. 4.1 Rigid Spheroid Model 

According to eq. (10) all the information a"bout the 
rotating nucleus is contained in I, Hence, most of the pheno- 
menological models aim at obtaining an analytical expression 
for I in terms of the shape parameters of the nucleus. Thus, 
for example, assuming the nucleus to be a rigid body of 
spheroidal shape, the value of I about an axis perpendicular 
to the symmetry axis is given by 

1r = (2/5) Mr 2 (1+ 0.31p + 0.42p^ + •••) (11) 

1 /5 

vdiere M is the mass of the nucleus and = 1 ,24 fm. , is 
the radius of equal volume sphere. It is found that the value 
of Ij^ is about twice that of 2° "to zero 

as p goes to zero. Thus, the nucleus cannot be treated as a 
rigid body, 

11. 4.2 Irrotational Flow Model 

Various other models have been developed assuming the 
nucleus to be a non-rigid body. The first of this kind of 
models is the irrotational flow model of Bohr and Mottelson 
(1953) based on the principles of hydrodynamics. They treated 
the vibrational-rotational mention of a deformed non-rigid 
liquid drop assuming the parameters p and V to be time 
dependent, and the flow to be divergenceless and irrotational. 
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The value of moment of inertia thus obtained for a general 
non-rigid asymmetric rotator with quadrupole deformation is 

Sin^ (Y + 27ro/3) (12) 

whgre 3 = 1 ,*2,3 correspond to the three principal axes and 
B is the mass parameter given by 

B = (3/8ii) (15) 

Bor axially symmetric (/= O) spheroidal nuclei (a, = a^ < a,) 

1 2 3 

we have 


"HI 


= I. 


H2 


= Ijj = 


and = 0 


(U) 


^H3 “ ^ means that the spheroid does not rotate about its 
symmetry axis. Although Ijj goes to zero as p goes to zero, 
the value of Ijj is only one- third of Thus, we have 




(15) 


and hence it follows that the nuclear flow is neither irrotational 
nor rigid . 

11,4.3 Centrifugal Stretching Models 

Following the above conclusion a number of rotational 
flow models have been proposed of which the centrifugal stretching 
model (Diamond et, al. 1964; Sood 1968; Gupta 1969; Trainer 
and G-upta 1971) is the most successful one. 
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The hasic idea underlying the centrifugal stretching model 
is that a rotating non-rigid body undergoes a centrifugal 
stretching, and the displacement from its equilibrium position 
is controlled by restoring forces. The idea is based on the 
experimental evidence that for low-lying rotational levels the 
moment of inertia I increases with angular momentum 1. 

II. 4.3a Model of Diamond et. al. (1964) 

Diamond et. al* (1S64) proposed a semi- das si cal model 
based on the assumption that the deformation (p) increases with 
the angular momentum (L) , thus leading to an increase in the 
moment of inertia (I) with I, In the final calculations they 
used the following energy expression 

E(I) = (•h2/2Ij.(p))I(h-i) + (0^/2) (Pj. - (16) 

along with the equilibrium condition = 0 to obtain 

the value of the superfluous parameter With this model a 

I 

good fit was obtained for rotational states in strongly deformed 
nuclei, assuming the relation lur p^ given by the hydrod 3 mamic 
model to be valid. However, the rotational states in the 
transition nuclei could not be fitted by this method with reasonabl* 
accuracy. Subsequently, various other forms of functions for 
the rotational moment of inertia I(p) have been used in eq. (16) ^ 
for obtaining bettor fits to the experimental data. Such | 

calculations indicated that the increase in p is not large , 

enough to explain the deviation from the L(L+1 ) rule. 
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11,4.513 Model of Sood 

Sood (1968) applied the concept of centrifugal stretching 
to the study of rotational-rihrational motion of a strongly 
deformed even~even nucleus, assuming the nucleus to be a simple 
classical rotator (point~mass) and the restoring forces to be 
harmonic type. This simple model was successful in predicting 
the energies of low-lying rotational states only for the nuclei 
in the middle of the region 150 < A < 190, 

11.4.3c Model of G-upta 

Gupta (1969) further developed the concept by considering 
an axially symmetric prolate deformed nucleus to consist of two 
parts obtained by dividing the nucleus by a plane containing 
the axis of symmetry. The center of mass for each of the two 
parts (see Figure IY.1) can be defined and the total separation 
r^ of the two centers of masses for equilibrium position is 
given by 

= (3/4) (17) 

where a^ is the semi-mao or axis corresponding to the ground 
state. The displacement (rj^ - r^^) from the equilibrium position 
due to centrifugal stretching was assumed to be balanced by 
harmonic restoring forces. The total energy B(Ii) of rotation 
which is a sum of kinetic and potential energies is given by 

E(Ii) = (-^^ 213 ;,) (Kl* +1) + (^g/2) " ^o^^ 
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■vdaere C_ is tiie force constant and I-r is t3ae moment of 

g ii 

inertia in the state of angular momentum L. To carry out 
the final calculations, the rotational energy obtained in the 
form of an assnnptotic series with I and 0„ as parameters 

O g 

was used. The model gave better fits to the low-lying rotational 
levels when compared to other models. Also, an attempt was made 
to calculate the energies of p -vibrational levels. However, the 
agreement with the experiment was not satisfactory, 

A quantum mechanical treatment of the model is expected 
to yield Improved results. This is done in Chapter lY and it 
constitutes a part of the v/ork presented in this thesis. 

The G-ovemor Model 

The semi- class! cal models of Soofl (1968) and Gupta (1969) 
give a non-vanishing value of I for spherical nuclei, -vdiich 
is contrary to the experimental observation. In order to include 
this. Trainer and Gupta (1971) modified the centrifugal stretching 
model of Gupta (1969) assuming that the effective mass contributing 
to the rotational moment of inertia lies outside a central i 

geometric spherical core which is rotationally' invariant. They 
refer to this model as the ’governor model’ , Assuming the | 

rotationally invariant core to be a sphere of radius equal to : 

I 

the semi-minor axis of a spheroidal nucleus (a^ = a 2 < a^) , 

the following expression for the rotational moment of inertia | 

Iq. was obtained (Pralcash 1973)? 
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= (V5) Mg (a^ + + 2a^) 


(19) 


wliere M is the effective mass participating in the rotations 

3 

and is given by 


— (,A-Tt/ O) O--] ( ) ■^0ff 


( 20 ) 


T/diere feff effective density in the region 

the rotationally invariant core, Notice that % 


for a 


1 


I j ♦ 


outside 
es to zero 


Both energy level and quadrupole moment data were used 
to obtain the values of the parameters involved in the final 
calculations. Significant improvement in the fitting was 
obtained for neutron deficient and 'transition nuclei. 


A quantum mechanical treatment of the governor model 
with a variable rotationally invariant core is expected to 
improve the results obtained by Trainor and G-upta (1971). A 
brief disci^osion to this effect is included in the last section 
of Chapter IV. 


II*4-,5 Variable Moment of Inertia (VMI) Model 

In viow of the uncertainty of the dependence of I^^ on 
, Mariscotti et. al. (1969) have developed a variable moment 
of inertia (Vm) model by treating the moment of inertia itself 
as a va,riable in the semi-classical expression 


E(I) = 


(«2/2Ij^) l(I+ 1) + (C^2) (Ij, - 


( 21 ) 
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where is the moment of inertia in the state of angular 

momentum 1 obtained by setting (9Ej/0Ij^) = 0. Since no 
model of moment of inertia is made, only empirical relation of 
the quadrupole moment with the moment of inertia is obtained. 

The VMI model has been most successful in predicting the 
energies of rotational states in the ground as well as in the 
excited rotational bands of deformed nuclei in the rare-earth 
and actinide regions. 

In order to fit the rotational spectra of the nuclei in the 
transitional region characterized by the limit 1,82 < E(4)/E(2)< 
2*23 within the framework of YML, it became necessary to use 
negative values of the parameters which led to the concept 

of ’phase transition' in nuclei (G-oldhaber and Qoldhaber 1970). 

The nuclei in this region arc characterized by considerably 
large E(2) values with a drastic reduction in level spacings for 
the hi^cr states in the band. The qualitative explanation 
given for the negative values of was based on the assumption 

that in the ground state these nuclei are not only rigid but also 
'brittle' , and hence these nuclei, when excited, first undergo 
some kind of rearrangement or 'melting' so that further excitations 
cost much less energy. The apparent phase change between the 
ground state and the higher states was qualitatively interpreted 
as a transition from 'super conducting’ to 'normal' state 
associated with the symmetry breaking of pair excitations. 





Recently, several attempts have been made to establish 
the equivalence of VMI model with other models, and to give a 
theoretical foundation to the VMI model. lEtie equivalence of 
the VMI model with the Harris formula has been shown by Klein 
et. al. (1970), An application of Vlil model for the study of 
of quasi-rotational spectra, and a derivation of the model -vdiich 
is general enough to account for its successful application to 
both deformed and spherical nuclei has been given by Das et. al. 
(1970a, 1970b and 1971). 
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C H A P T' E R - III 
DIEAMIG3 OP RIGID BODIES 

111.1*1 Definition of a Rigid Body 

A rigid body is a system of particles bound together by 
internal forces 'which always keep the distance between any t'wo •• 
particles a constant and act along the line joining Idle particles-, 
irrespective of the external forces acting on the body. 

Alternatively a system of n-particles subject to the 
constraints 


I _ Y~7Z „ A 2 , 7^ „ n 2 = constant (1) 


3 ' "'1 *"0 

r^3 n^ 


,th 


^i^i^i rectangular coordinates of the i particle 


■vdiere 
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and is the magnitude of the vector along the line joining 

the i and j particle, constitute a rigid body, (The rigid 
body is an idealization of the situation that occurs in practice 
and is useful when discussing the motion of a solid as a vdiole, ) 

111.1.2 Degrees of freedom 

A system of n independent particles will have 3n 
degrees of freedom. Since, by definition the particles in a 
rigid body do not move relative to each other, the most general 
motion that the body can perform is a rotation (3 degrees of 
freedom) and a translation (3 degrees of freedom), thus giving 
a total of six independent degrees of freedom for a rigid body 
and hence (3n-6) constraints. 

111. 2 AXIS OP ROTATIOIT MD AIT&ULAR VELOCITY 

Let the particles situated at O' , A and B satisfy the 
rigidity condition. Fix any one of them (say O') and rotate 
the system about this point through an angle d^ as shown in 
Figure III,1, 

The rigidity condition implies that the axis of rotation 
is the same for the particles at A and B and for any other 
particle located in the body. Thus, a rigid body moving with one 
point (O') fixed has an unique axis of rotation and an unique 
angular velocity w = (df/dt) about the fixed point. Further, 
it is easy to show that the angular velocity w of the rigid 
body is independent of the choice of the point fixed in the body. 



Pig, III.1 III. 2 

Rotation about a Pixed Point. Motion of a Point (P) in 

Space and Body Pixed Co- 
ordinate Syctems, 



Pig, III, 3; * Siler Angles' 



Also, dr = X r. 

111,5,1 Ooordinate Systems and Euler Angles 

The configuration of a rigid hody is always specified 
with respect to a convenient set of coordinate axes. The choice 
of coordinate axes can as well he arbitrary. But, usually, a 
proper choice of coordinate system will simplify the problem 
tremendously. G-enerally, the following systems of coordinate 
axes are useful in the study of the dynamics of a rigid body, 

(1) A system of coordinate axes fixed in space called the 

space fixed system, denoted by S(ZfZ) or simply S. 

(2) A system of coordinate axes fixed in the body called the 

body fixed system, denoted by B(x^X2X^) or simply B. 

(3) A system of coordinate axes chosen to coincide with 
B(x^X 2X^) at any instant of time called the instantaneous 
system of axes, denoted by B*(xJjx|x^) or simply B' . 

In Figure III, 2, 0 and 0* denote the origins of S and 
B respectively, P is any point in the body. Then 

Z = R + r . 

Taking the time derivative of this we have 

u = Y + w X r (2) 

where u and Y are respectively the velocities of P and Ot 
with respect to S and (w x r) is the velocity of P with 
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respect to B, 

Since w does not depend on the choice of 0^ in the 
body, 0' can as well be the centre of mass of the rigid body 

in which case, f is the velocity of centre of mass with 

respect to 0, further, if 0 and 0’ are made to coincide 
with each other and v/ith the centre of mass of the body then 
u is simply w x r. 

She various coordinate systems are related by certain 
spatial transformations. In discussing any spatial transfor- 
mation there are two distinct points of view. In the first, the 
transformation is brought by a coordinate transformation, while 

each point P of space, along with all physical entities asso- 

ciated with it, remains fixed. In the second, the axes remains 
fixed while the physical system is itself transformed, Ihrou^out 
this vrork the first of these viewpoints j;s adopted. 

In the case where the origins 0 and 0* are made 

to coincide with the centre of mass of the; body, it is always 
possible to express a vectorial quantity either in terms of 
the coordinates in S(xrz) or BCx^XgX^). The two are related 
by a coordinate transformation involving a rotation of any one 
of the- set of axes, for Ti\hich three parair'eters are essential. 

Ihe mosj useful of them are the Euler angles (<j), 0, f ) as 

A 

shown in "Eigure III, 5. 
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({) iiS the angle of rotation about the Z-axis which takes X 
to x]| and X to x^ and ^ is the corresponding angular 
•velocity along the Z-axis, 


0 is the angle of rotation about x^ axis which takes Z 
to x^ and x^ to x'^ and § is the corresponding angular 
velocity along the xjj-axis. 

Ip' is the angle of rotation about x^-axis which takes 
x^ • to x^ and Xg to and is the corresponding 

angular velocity along the x^-axis. 


III.3«2 Resolutions of Angular 'Velocity w in Various 
Coordinate Systems 

The angular velocity w has components ® ^ 


along Z, x]j and x^ axes respectively which form a non- 
orthogonal system oi axcs. It is possible to obtain explicitly 


the components (w^/Wgj w^) and Wy, w^) of w along 

the two orthpgonal system of axes b(x^3C2Xj) and s(XXZ) 

• * 

respectively, in terms of Prom Pigure III, 3, we get 





From these we get 


/w,l 


I 

-2 

= R 


\-3 




where R = P Q 


-1 


(5c) 


The transformation matrix R explicitly given hy 

t 

Cos llr Cos (|) - Cos 6 Sin (() Sin l| Cos Sin (|)+Cos 0 Cos Sin if 
R =|-SinTj' Cos <() - Cos 0 Sin ({) Cos T|r -Sin if" Sin (j)+Cos 0 Cos <j) Cos <{) 
Sin 0 Sin (|) -Sin 9 Cos (j) 


Sin l|r Sin 0 
Cos If Sin 0 
Cos 0 


(3a) 


The transformation relation (3c) is true for any arbitrary 
vector A (Goldstein, 1950). Hence, R is the general trans- 
formation matrix vihich transforms the components of a given 
vector A in S to B and vice-versa* It is easy to show 

•-jn 

that R is a. unitary real orthogonal matrix since R R = I 
and the transpose of R is equal to its adjoint. 


III *4, 1 Angula r Mo mentum of a Rigid Body in Motion 

If m. is the mass of i"^^ particle, r. is its 
position vector and u^ is the corresponding velocity, then 
the angular momentum (I) of a rigid body undergoing a rotational 
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motion about its centre of mass is 

- ^ - 

L = S 3C u^) 

i=1 

Using eq. (2) we have 

II = S m^[Jr? w - (r^.w)r^ (4) 

As with any vector equation, one may take the components of L 
in any desired system of coordinates, OSaus, for instance, 
in S(xrz) is given hy 


I 






( 5 ) 


with similar expressions for hj and In BCx^XgX^) 

system we have 


^1 = m^ I UI2 + ^13) - ^i1^i2^2 " ^i 1 ^i 3^3 ^ 

with similar expressions for L 2 and Ii^. 


III. 4. '2 Kinetic Energy of a Rigid Body in Motion 

Consider a rigid hody in motion. If m^ is the mass 
of the i"^^ particle and u^^ is its velocity in S, then the 
kinetic energy of the system is given by 


Using eq, (2) we get 

T = £ i m.(Y+^j_)^ = 2 “ T*(^^) + 2 ^ mj_. 

i=1 2 ^ 
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i.e*, T = i + (Yxw) ^ 2 + 2 “ w^r| - (w.Tj^)^ ]]. 

Clioosing the centre of mass of the hody to he the origin of B 
and using eq. (4) 


T 


1 MY^ + i L.w 

2 2 


m +0] 

trans rot 


(7) 


where M is the mass of the hody and the vectors 1 and w 
can he expressed either in terms of the components in S(22’Z) 
or B(x^X 2 X^), 

Thus, the kinetic energy of the system consists of two 
parts; (i) > the translation kinetic energy of the centre 

of mass in S and (ii) T^^^ , the rotational kinetic energy 
of the hody about the centre of mass. If one is interested 
in the study of only the rotational motion of the rigid hody, 
it is useful to choose the centre of mass as the common origin 
of hoth S and B in which case. 


® “ ^rot “ ^ 2 ''' ^2^2 ^3^3^ 

...( 8 ) 

where 1^, 1 * 2 ;) or (I^ , L 2 , B^) are respectively given hy 

equations of the type (5) or (6). Using the latter, we have, 


T 


1 


3 
Z 

2 l,k= 


"Ik 


where the inertia coefficients ^Ifc 
tensor I given hy (G-oldstein, 1950) 


'*'1 (3) 

are the elements of inertia 
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i.e*, T = -| MV^ + (Yxw) ^ E ^ ^ 

Choosing the centre of mass of the body to be the origin of B 
and using eq. (4) 






trans 


+ T 


rot 


( 7 ) 


where M is the mass of the body and the vectors L and w 
can be expressed either in terms of the components in sCZfZ) 
or B(x^X 2 X^). 

Thus, the kinetic energy of the system consists of two 
parts; (i) ^-t^^ans ’ translation kinetic energy of the centre 
of mass in S and (ii) » the rotational kinetic energy 

of the body about the centre of mass. If one is interested 
in the study of only the rotational motion of the rigid body, 
it is useful to choose the centre of mass as the common origin 
of both S and B in which case, 


T = T, 


rot 


1 

2 






■W^ 




1 




...( 8 ) 


where (L^;^ (^1 » ^2’ ^3^ respectively given by 

equations of the t 3 ^e (5) or (6), Using the latter, we have, 


T 


1 

2 


3 

2 

l,k=1 


where the inertia coefficients 
tensor I given by (Goldstein, 


^Ik ^1 ^k 

^Ik 
1950) 


(9) 


are the elements of inertia 



I 


( 10 ) 


S r^ul2 + xfj) 
“2 
-2 


■2 


-S 


2 -S “A 2^13 


■2 mj^x^2^^3 


2 m^Cxf^+xlg) 


In the above,, tne, diagonal elements are called principal moments 
of inertia and the off'-diagonal elements are called products of 
inertia, Ihus, there are only six indepatident inertial co- 
efficients and bhey depend both upon the geometry of the solid 
and the origin of thv^ '^rea. choaen^ and upon the orientation of 
these axBvS wT+h respect to the solid. Had we chosen to express 
the energy in terms of the coordinates in S(xrz ) , we would have 
landed up with equations exactly similar to (9) and (10) with 
(X,Y,Z) in place of (x^, X2» x^). 


In any solid, it is always possible to find a particular 
orientation of body axes, called the principal axes, for vdiich 
the inertia’ tensor is diagonal. In this case, 


^Ot = I ( + IgW^ + IjW® ) (11) 

w^ , Wg, w^ are given by eq. (3a). 

( Ij^ is written in place of I^^ for the case of principal axes.) 


111*4.5 Hamiltonian of the System 

In the absence of external forces, the rigid body 
undergoing a rotational motion has the Hamiltonian (eq. (8)), 



1 1 

...( 12 ) 

In paxticular, by choosing the B(x^X 2 X^) axes parallel to 
the principal axes of the system we get the most useful form of 
H given by 

H = 1 S (Ii!)A2It) (13) 

2 ^*=-1 3 D 

Inhere as given by eqs. (6), 

III,4.4 Conservation of L and H in case of the Motion 
of an Isolated Rigid Bo(3y 

We have shown that for a rigid bo<3y in motion with no 
external forces, the rotational kinetic energy is given by 
eq. (12) and the corresponding angular momentum components are 
given by eq* (5) or (6). In the absence of any external force, 
the system ¥)'ill not experience any torque and hence (dl/dt) = 0 
which means that I is a constant of motion. But the conser- 
vation of kinetic energy and individual components of I depend 
on the conservation of w and its components which, in turn, 
depend on the geometry of the solid. Here, we will discuss 
three interesting and useful geometries of the solid. We will 
use H as given by eq, (13). 

(a) Spherical lop t In this case we have I^ = I 2 = I^ = I. 

Hence L = I w which implies that w is a constant of motion 
and the Hamiltonian H given by eq, (13) takes the forn H = B /2I. 
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(t)) Symmetric Top ; In this case we have =12^= I^. Choose 
the Z-axis along the direction of I and the x-j-axis perpendi- 
cular to 1 so that the component 1^ of li along x^-axis is 
zero and hence = 0. So w must lie in the plane containing 
axes (x 2 »Z, x^) i.e,, in the plane of the paper as shown in 
Figure III, 4. 

If 0 is the angle between Z and x^-axes (refer to Fig. III,4)j 
then 

~ = If Cos 0 / I^ 

and W 2 = ( 12 / 1 ^ ) = Wp^ Sin 0=1 Sin 0/1^ 

hence, w^^ = (bAi). 

So w^ and are both constants of motion and hence w 

5 pr 

is again conserved. In this case eq, (13) fon H takes the 
form 

H = ( 1^/21^ ) + ( l 2 /2 ) ( 1/Ij - 1/1^ )■ (H) 

(c) Asymmetric Top ; In this case we have I^ ={= ^2 ^ ^3* 

^1* ('13) gives 

I^/I^ + + i^/Ij = 2E 

and = Ip-. 

Hence in B(x^X 2 X^), L moves along the path of intersection 
of the ellipsoid and the sphere defined by the above equations, 
which is necessarily a closed path. Hence, the motion of 1 in 
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B is periodic. It car. Be shown (landau and lifshitz, 1969) 
that 0 and "({r will have the same period, but ({) will be a 
combination of tv/o periodic motions. Thus, while 1 is a 
constant of motion, the top never exactly returns to its 
original position, 

111,4*5 Euler* s Equations of Motion 

It can be shown (Gioldstein, 1950) that the time derivative 
of an arbitrary 'vector in S is related to its time derivative 
in B, together with its vector product with 7/, the instanta- 
neous angular velocity of B, by the operator eq^uation 

(d/dt)g = (d/dt)g + w X (15) 

Operating (15) on the total angular momentum 1 , we have 

(dl/dt)g = (dl/dt)g + w X I) (16) 

Choosing B(x.|X 2 X^) along the principal axes of the solid 
we have 1. = I. w. and 

J J t] 

(dl/dt)^ = (di.j/dt)e^ + ((3l2/dt)e2 + (aLj/dt)e^ 

where (e^,e 2 ,e^) are unit vectors along the principal axes. 

In the absence of any external forces we have (dl/dt)g = 0 
and hence eq, (16) reduces to 

dw- 

, — + — Wt. w-] = 0 (17) 

dt I. ^ 

J 
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with, two more expressions obtained by the cyclic dombination 
of 0 ,k,l correspond to 1,2,3). Eqs. (17) give the 

description of the motion of an isolated rigid body with one 
point fixed in B and are called as Euler* s equations of 
motion of a rigid body. 


QUMTIZATIOI OE THE RIGIB BODY MOTION 
In III. 4. 3 j we have seen that the rotational motion of 
a rigid body can be described by an Hamiltonian of the type 
given either by eq. (12) or by eq, (13). Of these, the latter 
is relatively simple and convenient to handle. The solution 
of the associated Schrtfdinger equation will yield energy eigen- 
values and state functions which are very useful in the calcu- 
lation of other physical properties of the system. Hence, to 
study the quantum mechanics of a rigid body it is essential to 
know the total angular momentum (1) and its components in B 
and in S (1^, respectively) in their 

operator form and their various quantum mechanical properties. 


In III, 3.2, we obtained the transformation relations for 
trainsforming various components of w from a given set of axes 
to another. Thus, eqs. (3a) to (3b) imply 

• •• *■ 

(w.^ , Wg , W^) < 5>(({), 0, Ijf) « ^(w^;, W-y, Wg). 

111,5.1 Orthogonal Coordinates 

The Euler angles ({j),0,Tjr) are coordinates of the system 
referred to the set of axes (Z, x^ , x^^) which is non-orthogonal 
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(refer to figure III* 2), The SchrtJdinger equation for describing 
any system can "be written only in teims of an orthogonal (like 
the Cartesian) set of coordinates. In our study, both, 
and S(Xrz) form two independent orthogonal set of axes. If 
, a 2 » oc^ and ay, a 2 are the angles of rotation about 

the coordinate axes in B and in S respectively, these will 
then constitute two independent orthogonal set of generalised 
coordinates for describing the system* It is necessary to express 
these angles in terms of the Euler angles ((j)0ljr) for explicitly 
obtaining angular momentum operators which requires 

(a^ttgtt^) >(({) 0 f) ^(ayaya^) 

111,5.2 Angular Momentum Operators 

Since a^ya 2 ja 5 and ay, ay, represent angles of 
rotation about the coordinate axes in B and in S respectively, 
we have 

= a^ (d = 1>2,3) and Wy^y ^2 = “x,Y,Z 

and, according to eq. (11) 

1 ^ .2 

H = i 2 I. af 

2 3=1 3 D 

Thus, the quantities that are canonically conjugate to the 

a . * s are given by 
J 

aH/da_. = I. a. = I. w. = I. 

J J J J J J 
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The L. ’s are then represented in the Schrttdinger picture 

J 

by the differential operators 


L. = -i* (a/3a.) (18a) 

*J J 

Now, the partial derivatives with respect to a. are related 

J 

to those with respect to ({), 0 , ijj by the following equations 
( Appendix IE , 1 ) 

^ W * 'fef- - t oot 6 

5% = ili4 4“ ■ ’I’ ® ef" 

3-9 
3a^ - ^ * 

Using these in eq. (18a) we obtain Ii-j, I 125 ^3 explicitly 
in their operator form. 

Working along the saDie line with ajf we get 


"X,Y,Z 


“ (S/ScXj- -y g) 


(18b) 


Again, the partial derivatives with respect to 
related to those with respect to (|), 0 , I' by the follov/ing 
three equations (AppendixII.I ) 

^ = -Sin ({) co-fc e 4“ ^ ^ ■§1“ ^ W 

jL 

5 I- = Cos $ Oot e jI" + Sin 4) ^ - |jf-| 5 |- 

d ^ 8 

3a2 "W ' 


(19b) 
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Using "these in eq, (ISh), 7/e ge"ti explicitly in their 

operator form. 


Thus, we have for the operator 

I, = (I, .Lg.ij) = 

It is easy to verify that 


a 


a 


^aa^’day^aa^'^ 


£2 = + 1,2 ^ L^. 


j2 , y2 . t -2 


and each is equal to (Appendix H.l ) , 

““■■5“' ^ (Sin 0 + — 1 * ^ "I 

sin'^e aib"^ ® sin^ 0 at^ J 




. . . ( 20 ) 

Explicit expressions for (I-j ,l2,Ii^) and (Iy,Ly,I<2) terms 
of the Euler angles yield the following commutation relations 


-ih I^ and cyclic in 1,2,3 ; 
and [JlyjLy] = -fiS. and cyclic in X,Y,Z , 

The above are the usual commutation relations for angular 
momen"tum operators. It is easy to verify that the operators 
I , li 2 3 1^ ^ a»n. d 

commutation relations are listed in AppendixII, 2 , 


^X’^’^Z Hermetian, Other important 


111 * 5.3 Angular Momentum Eigen-Eimctions • 

In the last section we have obtained the components of 
angular momentum I in their operator form. In Appendix H, 2 , 
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we have listed tlie various conmutation relations satisfied by 
these operators. It is now essential to study the angular 
momentum eigen-functions and their properties in order to obtain 
eigenvalues and eigen- functions of the Hamiltonian given by 
eq. (13). 

From the commutation relations (Appendix II2) , it is 
clear that the opera;^ors ^3 commute with each 

other. Hence, it must be possible to find simultaneous eigen- 
functions |liMK> of these operators which render and 

I^ diagonal with eigenvalues 1(1+1 ) , M and K respectively. 
Thus , 


E^Ilmeo 

= L(I+1)?l^iLM> 

(22a) 

1^1 IMO 

= m. 

\UK> 

(22b) 


= Ml 

\mK> 

(22c) 


Our task, at present, is to find the explicit form of jLffiO, 

The Euler angles (<{),0jljl) are the dynamical variables 
which completely specify the orientation of the B-frame relative 
to the S- frame. The rotation of the B-frame relative to the 
S-frame can be generated by a rotational operator R(((),0,|f) 
which according to the definition of Euler angles is 

R ( (t),0,l{r) = R (f) R(©) R(<{)) 

The explicit form of the operator R has been obtained in 
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Appendix 11*5 (Elliott, 1958) 
R(^ >0 »^) == e 


-ie -r -it T 

‘z ^ IT ^ Z 

0 © 


where Ij.Iy ,^2 components of the angular momentum 

L with respect to the S-frame. Clearly, the operator R is 
unitary*' Since equation (25) contains only angular momentum 
operators 1 * 2 clear that the operator R((t),9,t) can 
only couple a i^ the usual angular momentum representa- 

tion to other 1^(IM*) with the same L value. 

Writing this formally 


Tjr (IM') = R((|),e,f) Tlr(LM) = s <lm' 1r('; ) 1im> f(hM» ) ( 24 ) 

M’ 

Here is the same function as f but of coordinates 

T'ef erred to the rotated frame (B-frame), -> stands for the 
three Euler angles 

Eoliowing Wigner, we define the matrix element on the 
right hand side of the above expression as 

ij,* (■»> ) = <IM' |R(t» )i™> (25) 

is the familiar rotation matrix. Hence, 

f (IM') = 4* (:?) tClM') 

Multiplying hoth sides by 4st(S> ) and using the unitary 
property of B^ma trices (since R is unitary). 
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Let R^,R 2»^3 Rj,Ry,R 2 be tbe Cartesian components of 

R and, R^, R^^ and R^,R_|^^ be the corresponding spherical 
components of R along B(x^,X 2 , 2 :^) axes and S(xyz) axes 
respectively. Now, consider a system of particles whose 
coordinates may be referred either to B or S, The rotation 
(a) of the system with the reference frame fixed is entirely 
equivalent to a rotation (-a) of the axes of reference frame 
during which the system does not change, provided the axis of 
rotation is the same in both the operations. Let L be the 
angular momentum operator associated with the rotations of the 
system of particles, and if be the wave function of the system 
referred to S, while is the same function but of coordinates 

referred to B, Since imposing the same rotation on the system 
of particles and on the axes of reference B has no effect on 
the description of the system, the infinitesimal change in the 
wave function of the system caused by such an operation must be 
zero (Elliott, 1957) > 

( R + L ) If^ = 0 (27) 

Also, since if refeired to S is independent of rotation R 
of B, we have 

R f = 0 (28) 


I.I.T. i ' f 
CENTRAL 
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Using the equation (26) in (28) we have 

0 = Rif = E R ) If (M') (29) 

XVllVX ^ 

R being the generator of the rotations of B, operates both 
on ( t) ) and on , Accordingly, since R is a linear 

differential operator 

0 = 5f = 2 [(R]4,(iI))+l4,(^)R]f^(IM') 
r T T 

= C5o)) 

(R )) denotes the result of operation of R on 

( >) ) and hence is a simple function of 3 ) and not an 
operator. 

We now study eq. (50) for particular spherical components 
of R •, either R^, R^^ referred to B or RQ,R_j_^ referred to 

S. 

First consider R^ (same as R^) where upon, because 
l|r_^ (L2I’ ) ijr^ (IM*), eq, (30) becomes 

■ l, 0 “') = 0 - 

Since ijl^ (IiM») are linearly independent, coefficients of 
ijr^ (LM* ) may be equated separately to zero. Thus 

Rote that the auxiliary function ijr^ do not appear in eq* (31a,), 
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Similarly taking the spherical components 

(related to the Cartesian components hy S-+-]= + ^ 

are used to obtain (Elliott, 1957) 

1 


T rCl-mM’) (M*+1 )-i 2 t 

^ ~ ^ [ 2 J ^,£‘+1 

In obtaining the above, we have made use of the .known eigenvalues 

■\ 

of the operators = (l^+il2) and the relation ^2 

(Rose, 1957). It is of interest to notice that the commutation 
relations of the operators R^^ are 


mmmt 




dlR- 


(33) 


or, in Cartesian form [JR^,R2]1= -iHR^ (1»2,3, cyclic). 

These differ from the customary relations by the appearance 
of a negative sign. Also, compare eq.. (33) v/ith eq. (21). To 
examine the. components R^, R^.| of R in S-frame, it is 
necessary to perform the sum over M* in the second term in 
eq. (30) explicitly because, the effect of L^, on 

ljr^(LM*) is not obvious. 

Doing this, we get 


The operator L operates only on the particle coordinates. 
Row, taking the spherical component Rq (same as R^) ir. (34) 

J, [(Rq 2it(V ).)f (M‘) ]-*K M t(™) = 0 
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and, expanding ijr(LM) again, we have 

^1 [ ) ) - MR. ) 

Equating coefficients of (LM^ ) to zero yields 


(M') 


= 0 


R 


^ ,(•? ) = -BM ) 


’0 


(31t) 


1 lowing the same procedure for R_^^ starting from eq. (34)> 
we get 


^+1 ^ 




[ 


(l:pM) (1+1+1 ) -jl/2 

1 




M+1 ,M 


(i>) 


(32h) 


■In obtaining the above expression, we have made use of the 
eigenvalues of the operator 1^ = (Ij+ily) and the relations 

^+1 = + Ta ^+' 

Erom eqs. (31a), (3111) j (32a) and (32b) it is easy to deduce 
the result of the operation of R^ on ). We find - 

as expected - that 


= Kl+Dh^ I^,(^ ) (35) 

f.omparison of eqs. (31a), (31b) and (35) writh eqs. (22) 
shows that the function ( "5^ ) is a simultaneous eigen- 

function of R^, R^ and R 2 and in the S-frame of reference 
it behaves like an eigenfunction of angular momentum 1, with 
Z-component M and in the rotated frame (B-frame), I^n ('^ ) 
is an eigenfunction of angular momentum 1 (but for the negative 
in (33)), with 3-component M‘ -vdiich corresponds to E in 
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|LM[l>. Important properties of 1^(5) ) are listed in Appendix 
II » 4" • 

In the rotational motion of a rigid "body, the operators 
H are associated with the rotational angular momentum of the 
body with components R^jRgjR^ in the B- frame and R^,Ry,R2 
in the S- frame. Thus, in the rotational motion of a rigid body, 
the angular momentum I associated with the motion of the system 
of particles is the same as the angular momentum R associated 
with a frame fixed in its body. 

In the case of a non-rigid body the total angular 
momentum associated with the rotational motion is a sum of the 
angular momentum associated with B-frame and the angiilar momentum 
associated with the intrinsic motion of the system of particles, 

111,5,4 Eigenvalues and Eigenfunctions of H 

Thus, having studied the angular momentum operators 
and their eigenfunctions, we are now in a position to solve 
the eigenvalue problem of a deformed rigid rotator quantum 
mechanically the Hamiltonian of which is given by eq. (13)* 

That is 

H =: 2 (l.V2I-i) (13) 

;)=1 

vshere I^ are the principal moments of inertia. 

We shall discuss the eigenvalue problem for two physically 
important cases of ellipsoidal ts^e of geometry. 
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(a) Synnnetric Top : The case of axially symmetric rotator 

is interesting to study as it is possible to obtain the solution 


in a closed form. Choosing the x^-axis in B(x^ ' as the 

symmetry axis, which implies = , the Schrb dinger 

equation with the Hamiltonian given by equation (15) is 
(Davidson, 1968) 


with 


E 

sym 


HlDMK> 

- 

- 


rKl'+D 1 

L T * ^ ^ 



(36) 

(37) 


Thus, the state functions of eqs, (22) are also the 

eigenfunctions of the Hamiltonian for the symmetric top. We 
have already shown that the state functions }IMK> are nothing 
but DjJ^(([), 0 ,f) . Various properties of these functions are 
discussed in Appendix 11.4. 


(" 6 ) Asymmetric Top (I .| =t= =}= 1 ^): A general rotator does 

not possess an axis of symmetry and hence K is not a good 
quantum number. Also, this is clear from the fact that althou^ 
commutes with and commute v/ith the 

Hamiltonian H for a general rotator given by eq, (13). 

However, since I 2 ; ^3 commute with each other, it is 

possible to represent the state functions |1M> of an asjTnmetric 
top as a linear combination of states of a symmetric top 
( Davidson , 1 9 68 ) , Thus , 


+I 

2 AT.ilM;> 
K=-l ^ 


IM> 


(58) 
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With 1LMK> as the basic set of wavefunctions , the eigenvalues 
of H can then he obtained by treating the asyiametry as a 
small perturbation over the symmetry which is a fairly good 
approximation for ellipsoidal type of solids. She Hamiltonian 
for such a system is 


H = 




t 2 

5 


1 


t 2 


+ (I/I 2 - 1/1^) (L^/2) = 


Tj Hh H 

sym per 


(39) 


The above expression is obtained by adding and subtracting 
the term in eq.. (13). Clearly for small deviations 

from axial symmetry I 2 ) the last term can be treated as 

a perturbation. Using |IMK> as the unperturbed state functions, 
the energy eigenvalues for an asymmetric top are given by 


E 



+ il^E 


vdiere B is given by 

AE = <mK\‘L^\mSi> ( I/ 2 I 2 - 1/21^ ) 

to a first order. The diagonal elements of 1? in |liMK> 
representation are 


= (^^2) (l(li+1) - £^) 


Thus, the energy eivenvalues for a slightly asymmetric top 
to a first order are 



E 




(40) 


Note that for I-| = I2J the above equation, gives the eigen- 
values for a synmetric top. Further details are discussed 
in the hook hy Davidson (1968), 

III, 5.5 Symmetry Properties of the State Function |Ui> 

For ellipsoidal type of geometry of the solid, the state 
function lliM> exhibit a number of useful symmetries. Most 
important of these is the rotational symmetry vjhich introduces 
arbitrariness in the choice of axes. By rotational 

symmetry we mean that under certain rotational operations all 
the physical quantities describing the system remain invariant. 
In fact, as we shall see, there are 24 such operations providing 
v;e use only right-handed system (there are also 24 left-handed 
systems). In order to investigate the symmetries arising out of 
such operations, we define the relabeling transformations , 

Tg and (Bohr, 1952); 

(i) produces the interchange 1 < — > -1, 3 < — > ~3. It is 
equivalent to a rotation of n about the 2— axis such that 



(ii) T2 produces the interchange 1 < — > 2, 2 < — > -1, It 
is equivalent to a rotation of 71/2 about the 3-axis 
such that ~ • 

(iii) produces the interchange 1 < — > 2 < — > 3 < — > 1 
such that = 1 . 
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Hence, "blie mos'fc general relabeling transforma ■felon tb.a'fe does 
not involve a change of handedness is T? tS 1?. Here p = 1,2; 

I 1 2 b 

q. = 1>2,3»45 r = 1,2,3; giving rise to 24 different principal 
axis systems « 

The symmetry properties of the state function jlM> 
under the operations T^ and T^ , ^ich have the eigenvalues 
+ 1 , are 

T^ lliM> = + |1M> 
lLM> == + lliM> . 

Thus, there are four classes of rigid rotator state functions. 
Since, the basic state functions in our case are the S 3 rametric 
top state functions, it is interesting to study the operation 
of T^ and T^ on |liME>. This leads to the following result 

( Appendix IE » 5 ) 

TJIM> Qin(Ii+K) jl,]^ 

It is useful to write equation (38) in the form 

lLM> = S fAir |1MK> + A T. in -K>1 (41) 

K=0 ^ -* 

Clearly the operation of TJ on |IM> shows that K can 

either be only even or only odd. Further, operating T^ 

on llM> and' relating the expansion coefficients in eq. (41) 
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with negative K to those with positive E vie get 
K - even 


(i) 


= (-1) Ag- 

state 

function 

is +|1M> 

(ii) 

^-K 

- -(-1)^ Aj- 

state 

function 

is -|1M> 

E - odd 





(iii) 


= -(-1)^Aj^ 

state 

function 

is -llM> 

(iv) 

^-E 


state 

function 

is +llM> 


From eq. (40) for the energy of an asymmetric top, it is clear 
that each level is characterised hy two quantum numbers 1 and 
K. Ihe above symmetry conditions on the state function (IjM> 
helps in obtaining the possible energy level structure in each 
of the above four cases. For example, consider the system 
belonging to case (i), Z-even class of symmetry. This system 
will have one level with L =0, no levels with L = 1 (since 
the state functions vanish), two levels with 1 = 2, one with 
1 = 3, three with 1 = 4 and so on. 
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CHAPTER - IV 

THE TWO MASS-PQIHT GEHTRIPUG-AL STRETGHIH& MODEL 

In Chapter II, it was pointed out that the two mass-point 
centrifugal stretching model of Gupta (1969) gives a satisfactory- 
description of hi^ spin states up to 1 = 14 in the gro-und state 
rotational hand of deformed even-even nuclei in terms of only- 
two parameters, viz., the distance (r^) between the two mass- 
points in the gro-und state, and the stiffness C of the potential. 
However, the treatment of Gupta (1969) was somev/hat semi-empirical 
and semi-classical, and the rotational energy was obtained in 
the form of an asymptotic series whose convergence was sensitive . 
to the magnitude of the total angular momentum (I), 

It will be shown in this Chapter that it is possible to 
obtain the energy of the rotational-vibrational motion of the 
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two mass—pointc in a closed form, by treating the two inass— point 
non-rigid system quantum mechanically v/ith suitable phenomenolo- 
gical potentials, Thus, use of a harmonic potential is found to 
yield an energy expression as a function of the parameters r^ 
and G, and is found to contain, in addition to terms occurring 
in eq. (180 of Chapter II, an L-dependent zero point vibrational 
energy term for a given nucleus. 

Essentially, the purpose of our present investigation 
is to examine whether the nuclear flow, which is known to be 
neither irrotational nor rigid, is ^equivalent' to the motion 
of two mass-points. 

The formulation, calculations and the discussion of 
the results are presented in the following sections. 

IT.1 EOBMUlATIOir 

The model is formulated by representing the rotational- 
vibrational motion of an axially symmetric prolate deformed even- 
even nucleus equivalently by the motion of two mass-points, 
situated at the centres of mass of the two fragments, with a 
phenomenological potential as a function of the distance (r) 
between them. The two fragments are obtained by dividiiig the 
deformed nucleus into two equal parts by a plane containing the 
axis of rotation and perpendicular to the symmeti^ axis. In 
such a representation, the increase in the moment of inertia 
of the nucleus due to centrifugal stretching resulting from an 



i 



f^rordal nucleus undergoing rotation and 
fugal . stretching ? 
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increass in "tliG roiaiional fnegnency is acccunied for, "by an 
increase in tlie distance betYreen the two' mass-points. Here, 
we consider only such Tihrations (p-t3rpe) Yiiich preserve the 
prolate shape and the volume of the nucleus ■under centrifugal 
stretching which means that the centres of mass of the two 
fragments will vihrate only along the nuclear symmetry axis, 
and hence no Y" -vibrations are considered. Such a system will 
have a reduced mass ( m ) equal to one quarter of the mass (M) 
of the nucleus and a moment of inertia (irp) given by (ha'vydov 
and Pillippov, 1957) 

Ij = llr^A (1) 

By simple geometry (Bigure ly.l ) it can be siiown that the total 
separation r of the two-centres of mass in their equilibrium 
position is given by 

where a^ is the semi-major aodLs of the spheroid in its 
equilibrium position and the subscript 1 is used to characterize 
the state of the system. Bor a prolate spheroid with quadrupole 
deformation we have (Bohr and Mottelson 1953) 

slj^ ~ Hq ( 1 y P Ij ) ( 3 ) 

Ihusj the equilibrium distance (^^j,) between the two mass— points 
is determined by the deformation and the total volume. Hotice 
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tlxat eq. (2) is independent of semi-minor axis 

For a rotating nucleus, ttie two-centres of mass would 
shift due to centrifugal stretching such that the amount of 
stretching is given by the displacement 3!he centrifagal 

force 

p = m hji = 1^/ (mr^) (4) 

causing the displacement (r^^-r^) are balanced, in our investi- 
gation, either by harmonic or by anharmonic restoring forces. 

In eq, (4) w is the rotational angular frequency of the nucleus. 

Alternatively, one can also consider the stretching of 
the nucleus in terms of individuail nucleons, ITaturally, when 
such a system rotates, due to centrifagal forces, the inter- 
nucleon distance would increase with the increase in w. 

If is the mass of each nucleon, the total centrifagal 

force F is 

A A 

P = Z m. r.v/^ = m w^ Z r. (5) 

i=1 ^ ^ i=1 ^ 

where again the integral over r^ defines the distance rj^ 
between the centres of mass and eq. (5) is same as eq, (4)* 

IY,2 CHOICE OP IHE POIMIIAIi 

A likely form of the potential is the curve A of 
Pigure 17. 2 which is characterized by a strong repulsion when the 
nucleus is squeezed in from its equilibrium position at r^ and 
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a po’ben’tial barrier "before fission. Aliliougli a more approprialse 
form of the potential to "be used is a single particle potential 
with two-centres of attraction each situated at the centre of 
mass of each of the fragments (Holzer et. al. 1969), we, for 
simplicity and in analogy with the model of G-upta ( 1969 ), have 
used a single particle potential with one centre of attraction 
as our interest, at present,, is only in the low— lying rotational 
states. iPor such states, we have approximated the actual 
potential in the neighbourhood of the equilibrium minimiua, 
by the harmonic potential (curve B of Bigureiy* 2) given by 

Vjj(r) = (s/2) (r - Tq)^ (6) 

■vdiere s is the stiffness of the potential. 

In order to investigate the ahharmonicity effect we 
have used the Morse potential (Morse 1929 ) which is extensively 
used in the study of diatomic molecules and is given by 

( 7 ) 

where the parameters D and r^ respectively correspond 
to the depth and the width of the potential (curve C of Pigure 
IV. 2), This potential has the general shape of curve A, though 
not the potential barrier feature, Hhile the right potential 
to be used for the entire spectrum is not the Morse potential, 
we have used it only to study the anharmonicity effect on rotational 
states wMch lie much lower in energy than the potential barrier. 




1 - 


exp \ - a(r - r’ 
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17.3 OHEORY AM ITS OOMPARISOU WITH OMR MOjDEDS 
17,3*1 - Witll Haxmonic Oscillator Rotential 

Using for example, the harmonic oscillator potential, 
the Hamiltonian (H) representing the relative motion of the 
system is 

H = (-hV2m) 92(r,e,<})) + (s/2) (S) 

■where r, 9, <[) are relative coordinates. The energy eigen- 
values and eigenfunctions of the system are obtained hy solving 
the SdarOdinger equation 

‘I’) “ B(1T,L) (r,e,<j)) . (9) 

where are the usual quantum nimahers characterizing the 

state of the system. Following the usual procedure, the total 
whve function is separated into the radial and angular 

parts by writing 

Since the potential is centrally symmetric, use of eqs, (8) 
and (10) in eq. (9) yields the following radial equation for 
the motion of the reduced mass (m) of the system 

vh.ere is the effective potential that the particle 

experiences and is given by 
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^y*3 THEORY MD ITS COMPARISON WIOH OWm MODELS 
IY,3,1 - With Harmonic Oscillator Potential 

Using for example, the harmonic oscillator potential, 
the Hamiltonian (H) representing the relative motion of the 
system is 

H = (-hVsm) ;7^(r,e,^) + (s/2) (ivrQ)^ (8) 

^lere r, 9, <j) are relative coordinates. The energy eigen- 
values and eigenfunctions of the system are obtained hy solving 
the SchrUdinger equation 


^•Ni T.M <!>) ~ E(lir,l) (r,e,<j)) (9) 


where R,I,M are the usual quantum numbers characterizing the 
state of the system. Following the usual procedure, the total 
wave function (T}r ^_^ ) is separated into the radial and angular 
parts by writing 


Since the potential is centrally ssomnetric, use of eqs, (8) 
and (10) in eq, (9) yields the following radial equation for 
the motion of the reduced mass (m) of the system 


[ 


2m 


dr" 


+ Vgff(r) - E(]J,L) J Ujj^Cr) = 0 (11) 


vdiere is the effective potential that the particle 

experiences and is given by 



V(r)(MeV) 



Fig. IV. 3 Plots of the harmonicC Vj-j) and the centrifugal 
terms and their sum (V^ff) as a function of the 
distance (r). 
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7eff(r) = . (12) 


A plot of the individual terms on the ri^^t hand side of 
eq. (12) end their s^an (7^^^) as a, function of the distance 
r is sho\m in Figure IV. 3. It is important to note that the 
effective potential is different for different I-values, with 
a characteristic equilibrium minimum at r = r-j^ for each 1. 


ilaking a Taylor expansion of ^^^^(r) aground the 


equilibrium position r^^ we have 


eff 


p(r) ^eff^^l^ 


8V 


eff 


8r 


1 9 

(r-rj^) + ~| (r-r3^)^+* 

r=rj^ ^ 2 


...( 15 ) 

Using the condition that — =0 at r = r and neglecting 
the effect of higher order terms we ha,i e 






Hence, eq. (11) becomes _ 
:2 


-7 

t 


% ^-(3q/2)(i^r^)2-E«; 


Ujjx,(r) = 0 


( 14 ) 


vdiere S' = E(H,1) - 7g^^(rj^) and = s + 


mr-T 


Eq. ( 14 ) is supplemented by the following boundary conditions; 
the eigenfunctions iijfj,(^) goes to zero for r -* 0 or «> . 
For lri» rj^ , sq. (14) approaches the Schx'ddinger equation 
for one dimensional harmonic oscillator and hence , the eigen- 
values of eq. (14) can approximately be v^ritten as 
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B' = (K + 1/2)11 f 


or (H + 1/2)11 V ^ ^ 3L(L+1 )il^/(mr^) ^CV^^) 


+ “ (rj-r^)"^ + 
2 ii o 


....(15) 


T/iiiere r^- = r + -“i — li-— 

1 o 3 

msr;^ 

and U, in general, is not an integer (Merzlacher 1961), 

However, for values of parameters characteristic of even-even 
nuclei in the rare-earth region, use of integral values of H 
leads to negligible errors. 

Eq. (15) gives, in general, the energy of vibrational- 
rotational motion of the two mass-point system, The groxmd 
state zero-point vibrational energy E(0,0) of the system 
is obtained by putting H = 0 and L = 0 in eq. (15). The 
energies of rotational states in the ground state rotational 
band is obtained by putting H = 0 in eq. (15) and subtracting 
out the zero-point vibrational energy E(0,0) from the resulting 
expression. The energy expression thus obtained is 


E(i,) = L(l+1) + I 


(16a) 


T/diere E^ = ('^)|^ ^ ( s/m)~ +" ~3Ii(I'+1 )h^/(m^r^) . - VT^mJ (16b) 


r^ + Ii(l+1 )fivCi’isr^) 


(16c) 



mr-? 

jj 


I 


m 

j. 
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(I6d) 


We ^low malce a comparison of tlie expressions obtained 
above with the corresponding egressions obtained ’ 037 ’ others 
on the basis of various models. 

A direct comparison of eq. ( 1 6a) with the energ 3 r 
ej:p)ression obtained b 3 r Gupta (1969) (eq. (18) of Chapter II) 
shows that 

(i) the two differ by the quantity E corresponding 
to the zero-point energy given by*^ eq, (16b) 

(ii) the energy E(l) is obtained in a closed form 
aided by the auxiliary eq, (16c) for the 
equilibrium position 

(iii) there is no brealcdown in rotational structure 
at high spin states. 

Farther, use of eqs. (2) and (3) in eq. (I6a) and a 
subsequent simplification, shows that eq. (I6a) is of the same 
fom as eq. (16) of Chapter II (Diamond et. al. 1964), with the 

difference that there is a zero-point energy term (E ) and 

2 

the moment of inertia I^, is nor/ given by mrj^ which is 
determined by values of the pai’ameters r^ and s. 

lastly, using r = (3/4)Rq (1 + )/ 5/4^ p) the 

instantaneous value of r (refer to eqs. (2) and (3) ) together 
with eq. (12) in eq. (II), the radial equation in the two mass- 
point centrifugal stretdhing model for the harmonic potential 
becomes 



74 





"i" 


7i' 


L(L+1) 



E(H,L) 


''ij\jjj(p) — 0 (17) 


T?4iere Br 




( 15 / 32 ) (3ml/8%) 

(45/128) (2i.IR2/5) (1 + 7 5/4^' p)^ 
s/m , 


(17a) 

(17b) 

(17c) 


in expression similar to eq, (17) has already been 
obtained, by Davydov and Chaban (I960) on the basis of the 
irrotational flo?; model of Bohr and Mottelson (1953) and is 


given by 



i 1 

dp^ 2Itj 

C -H 

1(1+1) + (p-p^)^ -eJ Uj^(P) = 0 

' (18) 

vdierc 

Bjj = 


(18a) 


% = 


(18b) 


%/% = 

s/m . 

(18c) 


A direct comparison of eq. (17) with eq, (18) reveals 
that oq, ( 17 ) is of the sane form as eq. (18), and that the 
specific differences in the values of various quantities may 
bo due to the difference in the specific flov/ pattom in the 
t?/o cases. Besides, it is rather interesting to note that 
although the mass parameter B^ is about half of B^, the 
two mass-point moment of inertia (1^) is about six times the 
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hydro djTiajnic value (Ijj)j and is about a third of the rigid 
spheroid value values of Also, notice 

that unlike 1^, does not go to zero with p. However, 

the frequency of vibration is the same in both the cases. 

Thus, the two mass-point model, based on the simple 
classical considerations of symmetry, represents a flow which 
is neither irrotational nor rigid , leading to the above values 
of the parameters. 


IY.3.2 With Morse Potential 

The expression for the eigenvalues with the Morse 
potential is well knovoi (Pekeris 1934, Pauling and Wilson 1935) 
and is given by 


(D-c^/ 2)2 2a(]>-c/2) 


E(N,L) = Cq - + 


where 


2^ Y2mrS+c^ 


^(N+ i) 


^0 = 




t 

aro 


+ 2 

t 2 t2 

aro a r^ 

_ § ) 

2 t2 ^ 


a r 


g 

2m 


1n2 


1-) 


...( 19 ) 


and 


= ■ ArC- 


gr 


t 


g r‘ 


ti 


21 


r- 1(1+1 ) 


0 -o 

^ T 2 

; where l’ = mr 


1, g, r' are the parameters of the potential. Bq, (19) 


differs from eq. (16a) for higher values of 1, The energies 
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of rotational states in the gromid state rotational hand is 
obtained hy putting IT = 0 in eq. (13) and then siibtracting 
out the zero-point vibrational energy E(0,0) obtained by 
putting IT = 0 and 1 = 0 in eq. (19). Ibe resulting equation 
is 


E(l) = 


°o - 


( l^-Ci/2)' 
1 +Cn 


a 


(I)-c^/2) a%i^ 


where S(0,0) = 


l/^aC5+c^ 
ahl 


2m 


- E(0,0) 


( 20 ) 




2^2 
a -XL 

Sm 


IY.4 OORBELAIIOIT OE THE IllRIITSIC QUADRIPOLE MOMl 

IHE MOIISITT OE IIIERTIA OIT IHE BASIS OE VARIOUS MOIELS 


A simple analysis of the available experimenta,! data on 
the basis of various models shows that the two mass-point model 
of Gupta ( 1969 ) successfully correlates the intrinsic quadrupole 
raoment da,ta with the moment of inertia data more satisfactorily 
than the rigid spheroid or the hydrodynamic models. 


The value of occurring in va.rious expressions is 

1 /3 

not well determined by a formula of the tyne R^ = 1*2 A ' fm, , 
since nuclei in the rare-earth region axe expected to have large 
skin depths (Hofstadter 1957) of the order of two femis. 

However> it is interesting to note that for spheroidal nuclei} 
the ratio of model moment of inertia (Ij^^) for the rigid spheroid, 
the hydro d 3 niami c and the two mass— point models > to the quadrupole 
moment given by eq, (6) of Chapter II is independent of R^, Thus, 



JABLE IY,1 




I { 


g... obtained for •various models usin^ experimental values of (l/Q) 


Nucleus 

"^expt 

(bams) 

^expt 

(10““^^ ke7 
sec^) 

Value of deformation obtained 
using experimental values 

of (I/a) 

Prp 


Ph 

1 52 

el®" 

5.85 

95.9 

0,372 

1.052 

2.003 


6.81 

142,5 

0.270 

0.728 

2,751 

■'54., 

6,08 

94.9 

0.390 

1.116 

1.924 

I56aa 

6,86 

131.3 

0.296 

0.806 

2.490 

^58sa 

7.30 

146.9 

0.282 

0,763 

2.626 


7.55 

155.1 

0.279 

0.754 

2.656 


6.17 

84.6 

0.464 

1.442 

1,668 


6.85 

118.0 

0.335 

0.925 

2.09 

^®l)y 

6.91 

134.7 

0.287 

0.778 

2.576 


7.13 

144.2 

0.278 

0.751 

2.666 


7.49 

159.2 

0,264 

0.712 

2.812 

68®" 

7.01 

115.6 

0.354 

0.988 

2,098 


7.23 

128.4 

0.324 

0.891 

2.278 

■'®®Er 

7.62 

144.9 

0.300 

0.816 

2.464 


7^.64 

146,4 

0.302 

0.821 

2*448 

'™I)r 

7.46 

147*9 

0.293 

0.795 

2.524 

■'®YD 

70 "“ 

7.39 

134.3 

0,312 

0.853 

2.367 

■''^Orb 

7.56 

138.7 

0.313 

0.855 

2.361 


7.77 

148.4 

0.301 

0,818 

2.456 


TaTjle 17,1 (Contd.) 



7.57 

152.7 

0.283 

0.768 

2.609 


7.40 

142.3 

0.307 

0.839 

2.401 

7?^ 

7 i 27 

128.5 

0.327 

0.901 

2.58 


7.37 

132.3 

0.325 

0.895 

2.270 

178 jjj 

6.78 

125.3 

0.317 

0.869 

2.327 

00 

t — 

6.73 

125.2 

0.319 

0.876 

2.311 

180^, 

74 

6,65 

114.5 

0.343 

0.951 

2.161 

182^ 

6,46 

116.7 

0.325 

0.895 

2.271 

00 

6,08 

105.4 

0.352 

0.892 

2.1 w 

186^ 

6.00 

95.4 

0.405 

1.175 

1,861 

76“® 

5.59 

85.1 

0.415 

1.215 

1.823 

1880. 

5.26 

75.4 

0.465 

1.450 

1.663 

'900s 

5.06 

62.6 

0.612 


1.367 
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^ere is a model dependent function, lor the tv/o mass- 

point model, Ij,T is given by eq. (17b) and corresponding 
model dependent function is 

% = /64)(M/Ze)(l + 0.63 P)V(P + 0.16 p2). (22) 

Expressions similar to eq. (22) can be obtained for the rigid 
spheroid and the hydrod 3 mamic models using corresponding model 
moments of inertia. However, such an analysis is not appropriate 
for the case of governor model since the effective nuclear density 
is treated as a parameter (Traiiior and Gupta 1971). 

Calculations are done for thirty- tv/o even-even nuclei 
in the rare-earth region for which experimental values 
(Mariscotti et. al. 1969) of quadrupole moments were available, 
IDieir moments of inertia are obtained from the Imown energies 
of the levels ?/ith 1 = 2. Using these vaJues of (I/Q^) in 
oq. (21), the values of p obtained for each model are tabulated 
in Table IV, 1. These values of p arc in turn used to calculate 
the values of for the two mass-point, the rigid spheroid 

and the hydrod 3 rnamic models, and are plotted as a function of 
the neutron number (A-Z) in Eiguro IV, 4. It is clear that 
only the two mass-point model gives satisfactory values of p 
and R^; For 90 < (A-Z) < 110, the values of R^ for the t?/o 
mass-point model are close to those given by the equation 
Rq = 1.2 fm. In this region, p c:' 0,35 and R^ Cif 


6.5 fln, , 
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for the two mass-point model, p 0.85 and ni 4*5 fm. , 
for the rigid spheroid model and p 2.2 and 2 fm. , 

for the hydrodynamic model. 

IF.5 OALQULATIOIS 
^^•5.1 Level Energies 

Por calculating the energies of rotational levels in 
the gromd state hand, we have used eqs. (16) in the case of 
hamonic potential with s and r^ as parameters, and eqs. (20) 
in the case of Morse potential with L, a and r^ as parameters. 

Por fitting the energy levels there were two alternatives. 
Either to adopt the least squares fitting procedure or to fit 
the low lying levels exactly to uniquely deteimine the parameters 
and then calculate the higher levels. While apparently, satis- 
factory fits to the levels of a large number of nuclei could he 
obtained with the least squares fitting procedure as has been 
done earlier by Mariscotti ot. al. (1969) and many others, the 
latter procedure is used to bring out systematic departures, if 
any. 

In actual calculations, in the case of harmonic potential 
eqs. (16) are used to fit exactly the lowest 2^ and 4*^ levels of 
about fifty even-even nuclei in the rare-earth region to uniquely 
deteimiine the parameters s and r^ for each nucleus, notice that 
is only an auxiliary parameter. These values of s and r^ 
are used to calculate the energies of levels up to L = 18. The 



maxim'uin departure of the calculated value from the experimental 
value (liiariscotti et. al. 1969) usually occurs for the hipest 
known level. Ihe values of maximmi percentage 

departure (AE)„ and the level to level comparison of calcu- 
lated energies are given in Table 17*2. 

Plots of variation of the stiffness C with the neutron 
number and the proton number are shown in Pigures 17,6 and 17,7 
respectively. 

The stiffness s’ characteristic of the equilibrium 
position r^ of the Morse potential is 

s’ = (a^ 7j^j(r)/ar2) = 2 D . (23) 

In order to stiidy the effect of anharmonicity, it is 
necessary to make a direct comparison of the results obtained 
using Morse potential with those obtained using harmonic potentia.1 
Such a comparison is possible provided v/e reduce the three inde- ■ 
pendent parameters D, a and r^ to essentially two, 

For low lying levels in the ground state the shape of the 
Morse potential closely resembles that of the harmonic potential 
and hence it is reasonable to assume that the stiffness s’ of 
the Morse potential to be same as that of the harmonic potential, 
thus reducing the number of parameters to essentially tv/o with 
the aid of oq. (23), The values of those two parameters are 
determined using eq, (20) and by fitting the known energies of 
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2*^ ajict 4 "^ s^taiies exactly. Tliese values of tlie parameijers 
are used "fco calculate tlie energies of rcfcafional levels up uo 
L = 18. aHie values of the parameters, the maximum percentage 
departure (AS) ‘ and a level to level comparison of calcu- 
lated energies with the experiment are given in faole IT,3» 
Although the ground state moment of inertia obtained with 

the Morse potential is about the same as that obtained with the 
harmonic potential, no noticeable improvement in tie level 
fitting is observed. 

Purther, the energy of the excited p— vibrational state 
(o'*** ) is calculated using the parameters needed to fit the 
lowest rotational levels exactly. A comparison of the calcula- 
ted vibrational energies with the eiQeriment is done in Pigure 
IY,5. However, a least squares fit to all the rotational levels 
and 0 level gives fits to within a maximum departure of 6 per 
cent as shown in Table If, 4. 


IY« 5.2 Quadrupole MoEient (Qqj^) and the Yalue of 
B(E 2 ; 4 “** - 2 ^)/B(E 2 | 2 '*' -* 0 “^) 

The reduced transition probability for transitions of 
the type L+2 - 1 is given by (Bohr and Mottelson 1953) 


B(E2; 1+2 -1) = 


g2 q2, 
° ‘^ol 


32r. 


^ 1+1 ) ( 1 + 2 ] 

21+3)121+5) 


(24) 


where is the intrinsic quadrupole moment of the nucleus 

in the state of angular momentum 1 given by eq. (6) of Chapter II 



84 

with the corresponding deformation for spheroidal shape. 

Sq. (24) can he used to obtain 

Bp . ? . ; . - g^) = (10/7 )(Qo4/Qo2)^ . (25) 

!Ehe value of the superfluous parameter rj^ obtained during 
the energy level fitting is used to obtain deformation 
using eqs. (2) and (3). A comparison of the calculated values 
of an^ "felie value of the ratio B(E2; 4"^ -* 2'^)/B(E2; 2*^ -* O"*") 

with those of the experimental values is done in lable IY,5. 

A comparison of the various results obtained above with 
experiments and other models is done in the following sections. 

IY.6 GOMPARISOE OE RESULTS mm EgEEOMlTS MB OSHER MODELS 
The experimental. values of the deformation p and the 
equal volume radius for a tj^ical spheroidal nucleus 

located around the middle of the region 150 < A < 190, are 
approximately equal to 0*3 and 6,5 fm, , respectively. It is 
clear from Table IY,1 and Figure 17,4 that the best values of 
p and Rq are obtained only for the two mass-point moment of 
inertia and hence the model correlates the quadrupole moment 
data with the moment of inertia data more satisfactorily than 
the rigid spheroid or the hydrodynamic models. 

In Table 17,2, for fifty even-even nuclei in the rare- 
earth region, we have tabulated the values of the parameters 
G and 1^, the experimental and the calculated energies up to 



L = 18, and the maximum percentage departure of the 

calculated energy from the experimental value for the case of 
hamonic potential, ’’^suallj?' occurs for the highest 

knoTOi level and is less than 1 per cent for 27 nuclei, less 
than 2 per cent for 7 nuclei and more than 2 per cent for 16 
nuclei, \diereas, for the 32 nuclei analysed by G-upta (I969)j 
( was less than 1 per cent for 17 nuclei, less than 

2 per cent for 7 nuclei and more than 2 per cent for 8 nuclei, 
considering levels up to L = 14 only. A least squares fit 
would give results as good as those obtained by ilariscotti et.al 
(1970), The values of moment of inertia obtained using the 
values of the parameter r^ in eq. (16d) are in agreement with 
the experimental values to within a maximum departure of 1 per 
cent, except for a few neutron deficient nuclei and **^*^0s nucleus 

The stiffness constants( C) that we obtain are about 
50 per cent larger than those obtained by Mariscotti et. al. 
(1970) and are considerably smaller than those obtained by Mosel 
and Greiner (1968). However, a number of prominent features 
characteristic to these nuclei as obtained by Mariscotti et. al. 
( 1970 ) on the ba,sis of the Variable Moment of Inertia model 
are also obseiwed for the two mass-point model. 

Ih*om Table IV, 2, it is clear that the values of the 
parameters and G vary rather smoothly from one isotope 

to the other. It is readily seen that, as a rule, large values 
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Calculated and Experimental Values of Energies, and tlie Values of 
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TiBLi! 17 « 4 

Values of the Harmonic and tlie Morse Potential Parameters and 

i S)xaax a least squares fit to all the rotational 

levels and E(o'*' ). 


Nucleus 

cni"^ 

( fm. ) 

(%) 

D 

(MeV) 

(10^^ cm“”* 

) (fin.) 

(%) 


0.39 

5x275 

4.6 

10i5 

1.4 

5.170 

5x0 

152sni 

0.38 

5.475 

5,9 

9.0 

1.5 

5,320 

6.1 

''5'^sm 

1.08 

6.455 

1.7 

45.5 

1.1 

6.380 

1.8 


0.39 

5.410 

4-7 

8.0 

1.6 

5.280 

6,2 ; 


0.96 

6.197 

2.6 

40.5 

1.1 

6.130 

2.7 


0.84 

5.886 

5.0 

36.5 

1.1 

5.S00 

4.0 


1.47 

6.175 

2.5 

52.5 

1.2 

6.110 

2.5 


1.46 

5.945 

2.4 

51.5 

1.2 

5.880 

2.4 ; 


2.02 

6.270 

2.4 

72.0 

1,2 

6.200 

2.0 ' 

1®Yb 

1.27 

6.015 

3.1 

53.5 

1.4 

5.950 

2.8 1 

■'™n) 

1.15 

6.028 

0.4 

48.0 

1.1 

5.960 

0.7 


1.11 

6.135 

0,6 

40.0 

1.2 

6,110 

0.5 ; 


1.81 

5.610 

0.5 

47.0 

1.4 

6. 1 60 

0.6 ’ 

1 

17% 

0.71 

5.750 

0.6 

36.0 

1 .0 

5.680 

0.6,' 

178Hf 

1 .50 

5.610 

0,8 

45,5 

1.3 

5.540 

0,9! 

182w 

1.50 

5.346 

0.5 

63.0 

1.1 

5,285 

1 

0.3; 

1®®0s 

1.21 

4.350 

4.7 

51.0 

1.1 

4.320 

5.o: 
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TABLE 17.5 

4 " 4 " / 4 " 

Coniparison of Calculated Values of B(E2, 4 *^2 )/B(E2;2 **0 

v/itli tile Eicperimental Values 


nucleus 

'^o2 

(■barns) 

B(B2; 4^-2’^)/3(E2^ 


Experimental Calculated 

Experimental* 

Calculated 


5.95 

5.22 

1. 11 

3.25 


6. 65 

11.07 

1.44 

1 .65 


6.15 

4.78 

1.41 

3.45 


6.91 

9.26 

1.42 

1,72 

158<}d 

7.20 

10.90 

1 .31 

1 .58 

''®I)y 

7.13 

9.08 

1 .43 

1 .64 


7.18 

10.05 

1.47 

1.56 

«°Er 

6.52 

3.00 

1.36 

4.17 


7.58 

9.75 

1.45 

1 .61 

''®Er 

7.60 

9.23 

1.50 

1.52 


7.81 

9.12 

1 .28 

1.56 

ISOjjf 

6.83 

3.77 

1.56 

1.66 

182^ 

6.40 

2.18 

1.42 

2.09 


*froiii ALou-Leila et al.(l97l)* Hucl. Phys. 1 T^ A, 663. 



Fig. IV. 5 Variation of' the calculated and experimental val 


of ECo**"') with the neutron number (A-Z) 




C ( MeV ) 
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of correspond to large values of C. Ho^vever, certain 

values appear to deviate from this rule. The glaring cases 
are and Relatively high values of 

Iq are ohserved for (A-Z) = 98 and (A-Z) = 108 isotones 
except for ^Hf, In particular, has the ma3:imum value 

for W and ^Os and decreases sharply for nuclei in their 
respective neighbourhood. The anomalous behaviour at (A-Z) = 98 
had previously been observed by Stephens et, al. (1964), -sdio 
speculated, as a possible explanation for this effect, that 
the pairing correlations are reduced (implying large moments 
of inertia) because of the large gap in the Hilsson diagram for 
such nuclei. More recently, Duckv/orth (1967) has shown that 
breaks are also seen at (A-Z) = 98 and (A-Z) = 108 in a 
plot of the double-neutron separation energy as a function of 
the neutron ntimber (A-Z)i 

Table 17,3 gives the calculated values of the parameters s, 
D, a, rl j the level energies and the value of ( AE) for the 
Morse potential, less than 1 per cent for 18 nuclei, 

less than 2 per cent for 8 nuclei and more than 2 per cent for 
23 nuclei. The values of the parameter are in general 

about 2 per cent less than those obtained with haimonic potential. 
The value of the depth parameter (D), although small, increases 
with the addition of neutrons and decreases with the addition of 
protons which is an expected behaviour. All other features exhi- 
bited by the Morse potential are very much the same as those of the 
harmonic potential. 
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03ius, no noticealDle improvement in the level fitting 

was ohserved "when the anharmonic Morse potential '.dth exactly 

same stiffness as that of the harmonic potential is used. For 

obvious reasons, much better fits are obtainable with the 

Morse potential when the parameters of the potential are treated 

as independent parameters obtainable by fitting the energies of 
“i" "{** 

2 , 4 and 6 states, 

A few striking facts may be gathered from the results 

given in Tables 17,2 and 17,3. 

(i) The largest value of stiffness (C) (but not 1^) 

1 RO 

occurs for Hf, an almost 'rigid rotor.' 

(ii) A sharp change in between 90 and 92 neutron 

numbers in Sm, Gd, Dj and Er isotopes results 
in a similar change in the value of the stiffness 
C. This is contrary to the conclusion drawn by 
Mariscotti et. al. (1970) on the basis of the 
7ariable Moment of Inertia model. 

The energy S(o'^ ) of the first excited ^-vibrational 
state calculated using the value of stiffness 0 needed to fit 
the low^est rotational states exactly Yirere on an average found 
to be lower than the experimental values by about 55 per cent. 

The variation of E(o’^ ) writh the neutron number (A-Z) as ' 

i 

shown in Figure 17,5 has more or less the same trend as the 
experimental value, except for the case of Er and Tb isotopes, 
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thus rendering a good support to the observed variation of G 
with the neutron number (refer to figure 17,6). It is also 
found that a least squares fit to the rotational levels and 
0 level gives results vdiich are in agreement mth experiments 
to within a maximum departure of 6 per cent as shown in fable IY*4. 

Table IT ,5 gives a comparison of the values of the 
intrinsic quadrupole moment of the first 2"*" state, and th^ ratio 
B(E2; 4"*" 2'^)/B(E2; 2"^ -* 0^) calculated using eq. (25) for 

about a dozen nuclei, with the e:cperimental values (Abou-Ieila 
and Baivrish 1971). Except for a few neutron deficient nuclei 
and the nuclei in the transition regions, the agreement appears 
to be quite satisfactory. 

IT, 7 DISOUSSIOI T 

1 1 ' . 7 » 1 T axiation of the Stiffness of the Harmonic Potential 
with the Feutron and the Proton Kjimbers 

The variation of stiffness C with the neutron number 
(A-Z) and the proton number (Z) shown in Figures IT, 6 and 
IT, 7 respectively, on the TAhole, show the same general trend 
as those obtained by Mariscotti et. al. (1970) and Mosel and 
G-reiner (1968). With the addition of neutrons, the stiffness 
increases sharply for neutron deficient nuclei and tends to 
decrease smoothly for neutron excess nuclei. This peculiar 
variation is rather interesting and may be understood in the 
following ways. 
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Since the p -vibrations, in our case, are of harmonic 
type, it is reasonable to expect ECo"^ ) to be proportional 
to the square root of G. A comparison of the points in 
Figure IY.5 v^rith the corresponding points in Figure IV, 6 for 
various nuclei appear to render a good support to the above 
argument. 

Also, comparing Figure IV, 5 with Figure 11,1, it is 
easy to notice that the variation of E(o'*'' ) with (A-Z) has 
a close resemblance with the variation of quadrupole moment 
with (A-Z), Using this, the variation of C with (A-Z) may 
be understood in the follov/ing way. 

Experimentally it is known that the intrinsic quadrupole 
moment increases with the addition of neutrons for the nuclei 
in the first half oi the region 150 K A K 190, and decreases in 
the latter ha,lf (see Figure 11,1 ) resulting in a corresponding 
increase or decrease in the deformation p and hence the 
moment of inertia nucleus increases or decreases. 

Whereas, a decrease in the value of moment of inertia (away 
from the rigid body value) makes the nucleus * softer* resulting 
in a decrease in the staff ness of the nucleus (neutron deficient 
nuclei) , the increase in the value of moment of inertia (approach- ! 
ing the rigid body value) makes the nucleus * harder* resulting ' 
in an increase in the stiffness of the nucleus, and vdien I ! 

approaches the rigid body value (neutron excess nuclei), G attain 
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the maxinrum value aud "becomes fairly constant thereon, Thus, 
the nucleus Hf , i^hich is lnio\m to "be one of the "best approxi- 
mations to an ideal ’rigid rotator’ has the largest value of 

1 80 

stiffness. The same result has "been obtained for Hf by 
Mariscotti et. al. (1970) and by Mosel and Greiner (1968) 
independently. 

The variation of the stiffness may also be understood 
in terms of the increase or decrease in the value of pairing 
energy (Mosel and Greiner 1968) when more and more neutron 
paifs are added to or taken out of the nucleus. 

The variation of stiffness C w^ith the proton number Z 
is shown in Pigure IY.7. Although, on the whole, the stiffness 
decreases v/ith the addition of protons, the trend for a given 
set of isotones, in general, is very similar to those obtained 
by Mariscotti et. al, (1970), the only exception being for 90 
and 92 neutron numbers. The stiffness is almost constant for 
(A-Z) = 90 and is approximately equal to 1,7 MeV. 

Such a variation could be attributed to the competition 
between the Coulomb and pairing forces in the nucleus when more 
and more proton pairs are added (Mosel and Greiner 1968), However, 
as we have not taken the effect of these forces into consideration, 
it is hard to draw any such conclusion. 
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IT. 7*2 Generalization of tlie Two Mass-Point Model 

We have seen that the simple two mass-point model 
exGplains most of the experimental data quite satisfactorily. 
However, the main drav/hack of the model is that the moment 
of inertia does not go to zero for spherical nuclei. In order 
to account for this, Trainor and Gupta (1971) proposed a 
governor model, in vdiich it is assumed that a central spherical 
por*tion of radius equal to the semi-minor axis of the spheroid 
does not participate inrrotations (rotationally invariant core: 
RIG) and obtained satisfactory results. It is rather interesting 
to note that the two mass-point model and the governDr model 
moments of inertia are obtainable as limiting cases of a 
generalised RIG model by considering the RIG to be of arbitrary 
radius (d) for a deformed nucleus of ellipsoidal ^ape 
(a^ > a 2 > a-j)« Ihe distance between the centres of mass of 
the two fragments, obtained by dividing the portion of the RIG 
nucleus outside the RIG into two equal parts, is given by 


.J _ 3 / ^1 ^"2 ^3 

^ 


a^ a2 — d 


5 


) 


( 26 ) 


vshere d ^ a^ < a 2 < Ihe corresponding generalised two 

mass-point moment of inertia is 




4 2 


-(27) 


where = M-M^ , is the mass taking part in the rotation, 
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Mg is the mass of the RIG and M is that of the nucleus. 

Clearly, eqs. (26) and (27) respectively go over to the 
case of the two -mass— point model and the governor model for 
d = 0, a^ = b . 2 \ and for d = a^ = a 2 < a^. Results obtained 
on the basis of the governor model (Trainor and Gupta 1971) 
are found to be in better agreement with the experiment* Hence, 
a quantum mechanical, treatment of the generalised two mass— point 
model is expected to yield improved results. However, we have 
not included this work here, as it is more or less a repetition 
of the case of the simple two mass-point model already discussed 
in detail. 
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CHAPTER - Y 

&EHBRAI. FEATURES OP COHTIHUUIjl MOEEH 

It is a widely accepted fact that the nuclear shell 
model provides a satisfactory explanation of a large number 
of experimental results for the ground and excited states of 
li^t nuclei (A < 40 ). But, for the nuclei in the region 
150 < A < 190 and A > 220 the shell model fails to account 
for the enormously large values of quadrupole moments, the 
properties of ground and excited states of doubly even nuclei, 
the magnetic moments of certain nuclei, the phenomenon of 
fission and so on. Such large effects may be expected to 
arise from the coordinated- both particle and surface- motion 
of a large number of nucleons inside the nucleus. Such a 
motion is usually termed as collective motion. 


!Eh.e close packing of the particles in the nucleus and 
the existence of a relatively sharp nuclear houndai^ have led 
to the comparison of a nucleus r/ith a liquid drop the energy 
of which is expressed as a sum of surface energy , volume energy 
and electrostatic energy (Uiels Bohr 1956? Niels Bohr and 
Kalckar 1957). Hius, an alternative to the shell model approach 
is the continuum model of the nucleus in which neither the nucleons 
themselves nor the interactions between them appear e:rplicitly, 

Ihe continuum model expjlains the collective phenomenon 
more satisfactorily than the shell model. 51ie model assumes 
the heavy nucleus to be a continuous) homogeneous) incoinpressible 
and uniformly charged liquid drop with a sha3?p boundary— a 
hydrod 3 mamic concept. The change in various properties of the 
nuclei can then be attributed to changes in shape and curvature 
of the charged liquid drop. Such an approach was first used 
by Niels Bohr and Iheeler (1959) to explain the phenomenon of 
fission. Subsequently, Rainwater (1950), Bohr and Mottelson 
(1952-55), Hill ond feeler (1955) and many others have used 
the hydrodynamic model for the study of deformed nuclei. 

The rotational flow model presented in this thesis is 
based on the continuum model of the nucleus. Hence, the 
general principles underlying the dsrnamics of a fluid drop 
and the a.pplication of these principles for the study of 
nuclear hydrodynamics will be reviewed in this chapter. 
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Y.1 DEFETIIIOII AIID BASIC PROPERTIES OF A gLUID 

A fluid may "be defined to iDe a continuous aggregation of 
molecules v/iiich. readily yield to the slightest stress applied. 

'Sie main difference "between a fluid and a solid is that a fluid, 
no matter how ■viscous, will yield in time permanently to the 
slightest stress, whereas a solid no matter how plastic, requires 
a certain magnitude of stress hefoi>e it will flow. A fluid may 
he either a liquid or a gas, the former is relatively incompre- 
ssihle vrhilst the latter is highly compressible. 

A fluid particle is an infinitely small volume element of 
the fluid Ti’ihich is large enough compared to the size of the 
molecule, hut small enough compared to the volume of the hody 
under consideration. Fluids considered in practice contain 
an infinite number of such volume elements and hence possess 
infinite degrees of freedom, 

A perfect fluid is a convenient idealisation of the 
normal fluid and is helpful in achieving mathema.tical simplicity 
of the equations involved in obtaining a solution. In many cases, 
it is fo-und to he a very good approximation to the actual fluid. 

A perfect fluid may he defined a.s one which is non-viscous, 
implying that there are no shearing stresses during the motion 
of the fluid, or incompressible, implying that the volume 
roma,ins constant, or both. 
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Y.2 BASIC lOmiAJjim 

The "basic pro"blsni in fluid dynamics is to find explicitly 
the velocity of a fluid particle in terms of its space-time 
coordinates. The formalism is based on the continuum model of 
fluid having appropriate continuum properties so defined as to 
ensure that on the macroscopic scale the "behaviour of the model 
duplicates the "behaviour of the actual fluid. In v/hat follows, 
we will discuss the laws governing the kinematics of fluid 
motion and the corresponding forces in fluids, and the dynamical 
equa,tions of motion. 

V,2,1 ’Lagrangian* and ^Eulerian* Modes of Description 

The motion of a fluid may be investigated by two 
different methods; the La.grangian method and the Eulerian or 
flux method. 

In the La,grangian method, we fix our attention upon 
an element of fliiid whose initial position vector (a) is 
specified at some given time t^, and follow its motion 
throughout its history. The coordina.tes x^ and t are 
themselves functions of the and t. Theii any quantity P 

associated with the motion of fluid will be a function of a 
and t. The velocity component v^ and acceleration components 
f^ are then 

v^ = 5x^/3t ; 


% “ 9"Vj_/at = (3^x^/3t^). 
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In the Bulerian method, v/e fix om? attention upon a 
particular point of space occupied hy the fluid and observe 
what goes on there. The velocity components v^ are taken 
as functions of and t, Then, any quantity i* associated 

V7ith the motion of the fluid is a function of two independent 
variables- the coordinates of a point in space and of time 

t. It is easy to see that the total change in P vdth respect 
to time t is 

3 

dP/dt = ap/at + z v.(ap/ax.). (1) 

i=1 

The d3?namical equations obtained by the Sulerian method 
are relatively simple to handle, particularly for the case of 
steady motion. To pass from the Euler ian to the lagrangian 
method, the equations (dxj^/dt) = v^^ , with initial conditions 
x^ = would have to be solved. Throughout this thesis the 
Eulerian method has been used, 

"^•2.2 Descriptions in the Body CB) and Space (S) Pixed Systems 
Por studying pure rotations, it is convenient to choose 
the centre of mass as the common origin of both B and S and 
the axes of B to be along the principal axes of the body under 
consideration. Purther, it can be shown (Goldstein 1950) that 
so. (15) of Chapter III is valid even if the body is non-rigid. 
The equation is 

(d/dt)g = (d/dt)^ + w X . (2) 
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Hie above eq_uation is of prime importance in obtaining tbe 
components of a rector in terms of tne coordinates either in B 
or in S* 

(a) Yelocity 

Operating eq, (2) on the position vector r of a fluid 
element v/e get 

u = v + v;xr (3) 

vdiere u and v are respectively the velocities in S and B. 

(b) Angular Momentum 

faking the vector product of eq. (3) v/ith r and inte- 
grating over the total volume of the fluid, v/e obtain 

I = i + R (4) 

vdiere L is the total angular momentum vector which consists 
of two parts, 1 , the intrinsic angular momentum due to internal 
flow, and R , the rotational angular momentum due to body- 
rotations, By internal flow v/e mean the motion of fluid elements 
v/ith respect to B, 

(c) Kinetic Energy 

The total kinetic energy E in S can be obtained by 
taking the scalar product of eq. (3) vvith u and integrating 
over the total volume of the fluid. Doing this we get 


E 


(5) 



i iy 


v&iere is the intrinsic energy of the system given hy 

= I 3^1 ^2 '^3 

and is the rigid body moment of inertia. 


( d) Acceleration and Force 

Ihe acceleration a in S of a fluid element can he 

s 

obtained hy the time differentiation of eq^, (3) and a subsequent 
mathematical simplification. Doing this, we get 


®'s “ % 2(w X v) + w X (w X r) (6) 

In the above, the last tyro terms on the ri^t hand side are 
respectively the Corioli’ s and the centrifugal terms. Also, 
eq. (6) gives the effective force acting on a fluid element of 
unit mass in the field of the remaining fluid elements. 

The set of eqs. (3-6) are most useful in the study of 
dynamics of a fluid body. 

V.5 PHYSICAD ASSUlgllONS MD BASIC BQUAIIOIS 

One of the fundamental physical assumptions is the 
conservation of mass leading to the ’equation of continuity’ 
given by 

(,9§/6t')+^»(fu) = 0 (7) 

yhere f is the fluid density, A fluid with an uniform density 
throughout its volume is said to be ’homogeneous*, Further, if 



tlie density’" is constant in time for an homogeneous fluid, then 
eq, (7) reduces to 
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V • n = ^ • V = 0 

the familiar divergence less condition, The principle of 
conservation of mass requires the constancy of volume of the 
fluid under consideration and such a fluid is said to he 
’incompressible’ . We will be dealing only v/ith such fluids. 

Jlnother important property of a fluid is the viscosity, 
property vdiich 8,rises basica.lly out of the friction between 
the fluid particles in motion and results in the dissipation of 
kinetic energy. The most general formulation of classical flu.id 
is completed by the inclusion of viscous effects. However, 
in our treatment, for simplicity, we have assumed the viscous 
effects to be negligibly small and hence, v/ill not discuss 
this property here. 

The third important property of a fluid is the ’vorticity’ 
or ’ circulation’ . It gives an indication of the rotational 
character of fluid in motion at a given point in the body of the 
fluid. 

The vorticity at any point inside a fluid with respect to 
S is the curl of the velocity field at that point and is denoted 

"by f 

% = 


"7 X u 


( 9 ) 



Thus, the vorticity represents for a deformable medium the 
generalization of the concept of aiigular Telocity of a rigid 
body and is sometimes called as 'fluid rotation’. Using 
eq, (3) in eq* (3) we hare 

^ = V 5: u = V V + 2w (10) 

% 

The ' flux of Torticity* through the area bounded by the fluid 
is defined as 

//f.dX = //(vxu)*dl 

vhere dA is an elementary surface area. 

■^.pplying Stoke’ s theorem to the aboTe 

ff \ • dX = § u*dX . (11) 

The quantity on the right hand side of eq. (11) is defined as 
the ’circulation’. Thus, the circulation along a closed path 
is equal to the flux of vorticity through the area bounded by 
the closed path. It can be shown (landau and lifshitz 1959) 
that in a perfect fluid the circulation around a closed path 
is constant in time. This is called the •la?/’ of Conserration 
of Vorticity’ , Thus, for a given velocity field, the vorticity 
is preserved for all time 3,t each and every point in the body 
of a perfect fluid. 

There are two t3rpes of motion of a fluid- the * steady* 
motion and ’unsteady or turbulent’ motion. The motion is said 



to iDe steady only vdien the velocity at each poin'G in the field 
is independent of time. Othermse, it is said tc he unsteady. 

In a steady motion the lines of motion coincide with 
the paths of the fluid elements and will never intersect each 
other* In majay problems the steady flow is a good approximation 
to the actual flow and helps to achieve mathematical simplicity 
of equations to be solved, 

lastly, since liquids always possess a ohai'p boundary, 
specification about the boundary will provide us with additional 
equations which are often useful in obtaining a so?rjition of the 
problem under consideration* 

^•5*1 Euler* s Equations of Motion 

Ihe equations of motion for a fluid arc obtained using 
Ifewton’ s second law of motion. Euler* s eijuatio.is of motions are 
obtained by equating the acceleration at any point inside the 
fluid to the sum total of forces acting on unit mass of the 
fluid situated at that point. !Ehe forces acting at any point 
inside the fluid can be classified into twoj the internal 
or the body force and , the external or impressed force, 

!rhe acceleration at that point is given by eq, ( 1 ) rath' u in 
place of E. Hence, Euler* s equations of motion are i 

(3u/3't) (r* 7 ) u = E^ + E^ • ( 12 ) 

In the above, given Ej^ + E^ , one can solve for n and vice-verE 



Y.4 lAJURB Qg FJCLEAP. kiUID 

It was pointed out in tiie beginning of this Chapter 
that deformed heavy nuclei exliihit strong collective properties 
v/hich are best explainable through a collective node! approach 
vherein the nucleus is approximated to a liquid drop with a 
sharp boundary and the energy of which is written as a sum of 
volume energy, surface energy and electrostatic energy. 

A consideration of the basic properties of nuclei wall 
lead us to make some meaningful assumptions about the nature of 
nuclear fluid. 

The nuclear density ( ]"' ) is known to be fairly constant , 
throughout the volume of a given nucleus except for a slight | 

tapering near the surfe,ce and hence, for a given nucleus, it is 
assumed that df/dt = 0 and tbe volume remains constant. 

This moans that the nuclear fluid is homogeneous and incompre- ' 
ssible . 

The basic problem is to obtain the explicit form of 
interaction between the nucleons inside the nucleus in the 
absence of external forces. The internal forces which make a i 

j 

significant contribution are those arising out of nucleon and ; 
electrostatic interactions between the nucleons and the protons I 
respectively, Purther, as the nucleus does not dissipate i 

I 

energy, the frictional force between the nuclear fluid particles - 
is assumed to be absent meaning thereby that the nuclear fluid is | 



’ inviscid. * Also, a study of the available quadrupole moment 
data (see 11,2,3) sbov;s that a deformed nucleus, in the region 
of interest, can be assumed to be either spheroidal or ellip- 
soidal in shape. 

Ihe above assumptions about the nature of nuclear 
fluid lead us to the following equations which form the basis 
for the study of nuclear hydrodynamics, 

d f /dt = 0 (homogeneity) 

V • u = V • V = 0 (incon^ressibility) 

= '9'xv + 2w = '^xu (vorticity) ( 13 ) 

6u/dt + (u*v)u = f. (absence of external forces) 

. (^j_/3-^) ~ "I (ellipsoidal boundary) 

In the above a^’s correspond to the principal semi- 
axes of the ellipsoid and correspond to the sum of effective 

nuclear and electrostatic forces at a given point in the body 
of nuclear fluid. 

Thus, we conclude that the nuclear fluid has the 
properties of a perfect fluid- homogeneous, incompressible and 
inviscid- with a sharp boundai^. Hence, the study of nuclear 
hydrodynamics should be relatively simple and strai^t forward, 

Y.5 (TYPES OP PLOW SMSKiLLLY GOHSIDEKEI) 

Although the nuclear force is known to be short-range 
strong and attractive, the explicit fom of interaction is not 
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yet knovfli. It is possible to obtain tlie effective force in tiie 
continuum model only when we knov/ the velocity field u 
explicitly. On the other hand, one can make some physicadly 
valid assumptions about the nature of nuclear flow to obtain 
the velocity field u explicitly and subsequently use it to 
obtain the kinetic energy of the continuum and the effective 
force, Ihe validity of the velocity field used is then confirmed 
by a direct comparison of the findings of the theory with 
experiment , 


The vorticity condition leaves us v/ith only two possible 
types of nuclear fluid flov/, the first of which is a flow for 
^'diich the vorticity ( ) is zero called the ’ irrotational or 
vortex free* flox? and the other for v/hich the vorticity is non- 
zero called the ’rotational or vortex* flow. Ihus^ 


I ^ X u 1 = 0 


* irrotational flow’ 


^ X u ~ ^ 


’rotational flov/’ 


(14) 


It is found that the nuclear moment of inertia obtained 
by assuming that the whole of the nucleus takes part in the 
rotation with a rigid flow (v = O), is about twice the experimental 
value (see II. 4.1) which means that the effective mass taking 
part in the actual rotation is only about half of the total 
mass of the nucleus. Hence, it is necessary to assume flows 
v/hich effectively decrease the value of model moment of inertia 
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from the rigid body value. Again, such flows can he either 
irrotational or rotational type. 

^•5.1 Irrotational glow 

The nuclear fluid flow T^hich is ir 2 ?otational ’with 
respect to S ^ould satisfy the follov/ing conditions 


V • n = • Y - 0 

"^xu = '^xv + 2w = 0 


(15) 


nearly u can he written as the gradient of some scalar 
potential function <{) called the velocity potential. Then, the 
velocity u is given hy 


u = - V <l> • ( 16 ) 

Comhining eq. (16) v/’ith eq, (15) we obtain 

l> = 0 (17) 


which is the familiar laplace equation for (ji. It is alv/ays 
possible to solve eq. (17) for a given boundary condition. For 
strongly deformed niiclei the ellipsoidal boundary condition 
may be used. Thus, knowing ^ , one can explicitly knov; the 
velocity field u and hence the energy of the drop and the 
nature of effective force. Hov/ever, the value of moment of 
inertia thus obtained (Bohr and Mottelson 1953) is found to be 
only one-fifth of the experimental value meaning that the effective 
mass taking part in the rotations due to irrotational flow is 



much less "bhaja the actual nuclear mass contributing to the 
rotations. 

Y.5.2 Rotational llovf 

We have seen that the value of nuclear moment of inertia 
lies between those predicted on the basis of the continuum 
model with irrotational flow and the rigid body model, and 
hence the nuclear flow is neither irrotational nor rigid, 
nevertheless, it is possible to postulate other rotational flows 
(rigid flow being excluded) in order to reproduce the observed 
values of moments of inertia of the nuclei considered. A 
detailed discussion of flows of this type is given in Chapter VI 
•vdierein we have proposed a rotational flovir with a constant 
vorticity for the study of rotational spectra of deformed even- 
even nuclei in the rare-earth region. 

V.6 SEMI-GLASSICAX TRjBAHMg 

Ihe usual treatment is that one begins with the classical 
theory of a continuum in order to obtain the total energy or the 
collective Hamiltonian of the system and a subsequent quantiza- 
tion of the system by usual methods. In such an approach, 
fundamentally, it is necessary to postulate the nature of fluid 
motion within the continuum of loiown shape. In an actual 
problem as a first step one attempts to obtain the functional 
form of the energy and angular momsitum at a given point inside 
the fluid drop by postulating a local velocity field at that point. 


iHie total energy and total angular momentum as a function of 
shape parameters of the drop are obtained by integrating the 
corresponding local quantities over the volume of the continuum, 

!Phe resultant Hamiltonian is vhat is termed as the ’ Collective 
Hamiltonian' and is subsequently quantized by employing the 
usual quantum mechanical rules, fhe Schro'dinger equation for 
the quantized Hamiltonian is then solved to obtain the eigen- 
values and eigenfunctions of the system. It is always conve- 
nient to quantize the rotational part of the collective Hamiiti. 
in the angular momentum representations, Oiie energy spectrum 
such a Hamiltonian depends on the inertia and the stiffness 
^ -^-^ettrs 5 and quite -different nuclear spectrum is obtained 
for different values of these parameters, "^n the study of 
nuclear rotations, the success of a model is determined mainly 
by its ability to give the proper functional fora of nuclear 
moment of inertia which can satisfactorily reproduce the experi- 
mental values. In fact, the extent of agreement of all theoretical 
predictions made on the basis of a model with experiments is a 
measure of the accuracy of the flow chosen and the quality of 
the Hamiltonian obtained. 
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GHAP2ER-VI 
ROPAglONiX HiOW MOPED (REd) 


It is well kno\m that the empirical value of the ground 
state effective moment of inertia calculated using the usual 
L(L+1 ) rule for the defomied even-even nuclei are given neither 
by the rigid spheroid model nor by the irrotational flow model 
of Bohr and Mottelson (1953). In fact, the representation of a 
spheroidal nucleus by a two mass-point system, with and without 
a rotationally invarient core (RIG) (Trainor and Gupta 1971 1 
Gupta 1969) are essentially attempts to parameterize the nuclear 
flow which, within the framework of the macroscopic model, is 
neither rigid nor irrotational. Krutov (1958) has used the 
concept of rotational flov/ by defining the collective rotational 
motion as the change of mass density distribution in time such 
that any motion ydiich does not change the mass density distribution 
does not contribute to the energy of collective motion. Rov/e (1970) 
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lias iii1;erpreted the assiomption the RIO as the existence 
of a rotatioiially invariant snpsrfluid core which remains 
as an inert spectator to the rotational flow of the remainder 
of the nucleus (the normal fluid) thus, referred to as the tv/o- 
fluid model. He has shovn that the adiabatic no del ^ which 
assumes the separate existence of rotational and intrinsic 
particle excitations, with very little coupling between them, 
implies rotational flow, vdiich is nothing but the rigid body 
flo ?7 except that the particles are not frozen in position but 
are free to execute any independent motion. Such a flow is 
found to lend support to the existence of clusters which 
participate in rotational flow as if they were elementaiy 
particles with their masses concentrated at their centres of 
mass thereby reducing the moment of inertia from its rigid 
body value. Recently, dupta (1572) has shoTOi that the linear 
dependence of the moment of inertia on the deformation 
parameter (p) is the characteristic feature of a system follow- 
ing the rotational flow motion and the RIG hypothesis, like 
the quadratic dependence is that of a h3rdro dynamical system. 

In view of the above and in viev/ of the fact that the 
nuclear flow is neither irrotational nor rigid, in this 
Chapter, we have attempted to propose a rotational flow model 
(RHI) based on the principles of continuum mechanics discussed 
in Chapter T. She model is developed using vortical flows 
in throe dimensions obtained by extending the irrotational 
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flow model of Bohr and Mottelson (1953) to ino^uce rotarioiial 
flow with finite but constant vorticity. fne rigid spheroid 
and the irrotational flov? models are obtained aj limiting 
cases of the proposed RBI, The model is fo^mJ to give more 
satisfactory results and is found to have many important physical 
implications. The outline of the model, theoretical fornula,tien 
and its physical significance, calculations anJ results are 
discussed in the follovang sections. 

TI.1 MOBEB 

The RBI is essentially based on the idea of extending 
the irrotational flov/ model of Bohr and Mottelson (1953) to 
include constant vorticity flows. 

We have seen in Chapter II that the energies of the 
low-lying rotational levels of deformed even-even nuclei in 
the rare-earth region are approximately given by 

E(l) = 1(L+1) (1) 

#iere B = 0, 2, 4, etc., and I is the nuclear i-otational 
moment of inertia in the state of angular momentur.' B. 

Ror a homogeneous, inviscid and incompressible fluid 
system with iDarotational flow we have y ‘U = 0 and ^x u = 0, 
where u is the velocity field in S. The corresponding moment 
of inertia (Ijj^) for a nucleus of ellipsoidal boimdary with 
semi-axes a^, a^, a^ (i, k = 1, 2, 3 and cyclic) with 
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ijii 

respect to tlie i principal axis is 


“■Hi 


= I 


Ri 


(a? - + a-j?) 


2 ) 


viiere a-, are given "by ea. (o) of Cli 2 ,puer II , ?jid 

^Ri rigid body moment of inerti,:, wdtli respect to the 

"t/ll 

i ■' principal axis. 

S:-:perimental data indicates that 


< Ir 


(5) 


aiiere ^gxpt obtained from ea. (1) using the energy of 
1=2 stats in the ground state band. 

It Y/as sho\m in CJhapter III that the velecit^r field 
(v) in B is related to u in S by the equation 


u = V + Yf X 


(4) 


Y<h.ere w is the angular velocity of B, Talcing the curl of 
oq. (4) v/e get 

Vxu = Vxv + 2w (5) 


Bor irrotational flow eq. (5) reduces to 

^x u = 0 ; and v - ^ = -2vr (6) 

for rigid flow, it leads to 

V X u = 2w j and ^xv = 0 (7) 



since y is zero. Eqs. (4), (7) and (8) su~r;e 2 t that it is 
possible to reproduce the values of , Itt end It, 

0Xpl/ ll s\ 

assuming a flow with a consta:it vorticity condition, namely, 

Vx u = - 2^ (6) 

Combining eq. (5) with eq, (8) we obtain 

X u = 2(w - (9) 

where Jv is a constant vector. ITote that o gives 

the rigid flow (Ig^) and w gives the irrotationa.1 flow 

(ljj)» We call ^ as the ’vorticity vector’ . 

Eor a system performing pure rotations, the potential 
energy of the system does not change, and hence the Hamiltonian 
is given by just the kinetic energy, Eie Hamiltonian H and 
the angular momentum I of the continuum with density f are 
given in S by 

H = ~ / f u^ ar ; L = / (r :c u) cr (10) 

2 

In the above, the integrations are over the entire volume of 
the system. By replacing u by v wc get similar expressions 
for the Hamiltonian and angular momentum 1 as measured 

in B, jQie relations between the quantities defined in S 
and in B, obtained from eq. (1) are 
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r- 




r (Lf - if)/2ip. 

A — ILJ- 


(i: 


i=1 


v^aere i = 1,2,3 correspond to tlie principal anc-s of the 
system; and w^. 

Essentially, the problem at present is to obtain a 
velocity field which satisfies the vortici'jy cosidition giver-i 
by eqs, (8) and (9) together v'ith the incompressibility 
condition, ^•v = ^»u = 0. Ihen the total rotational 
kinetic energy could be obtained using this velocity field 
in eqs. (10), (11) and (12). 

^1*2 IHBORBgICAIi POBMULAIIOR 

The constant vorticity condition and the divergencelcss 
condition in RH€ are analogous to the basic differential la;vs 
of magnetostatics. Hence, the most general form of the velocitjr 
field must be obtairiable, in principle, by solving the equations 
V X V = -2Ja; and .v = 0 for a specified boundary, Sie 
usual method of approach is to exploit V ‘V = 0; end solve the 
Poisson equation for the vector potential of v by G-reen*c 
function technique and hence obtain v explicitly in its 
analytical form (Jackson 1965), Howevor, we dispense with 
this method, as there are simpler ’ways of proposing the velocity 
field satisfying all the required conditions. 


^I*2*1 G-eneral Eoimulation in Ihree' Dimensions 

A rotational flow in three dimensions, in general, is 
characterised by the fact that each of the vectors v, u,5fi- and w 



vdll have non-vanishing components along the coordinate axes 
in B and in S, We look for a particular t;/pe of velocity 
field vdiich ■would satisfy all the flo'w requirements in general, 
and vdiich ’.vould reproduce the rigid hody flow and the irrota— 
tional flow under appropriate limits of the vorticity vector. 
One such flow in B satisfying eqs. (8) and (9)? and the diver- 
genceless condition for an ellipsoidal boundary is given by 


^i 


= 2a, 


to. 




/ (a| 


a^) —ft-- 2 




/ (a| 


a^) 


...(13) 


vdiere i, 3 , k= 1 , 2,55 with two more equations v/ith i, 3 , 1 : 
cyclic. In fact, eq* (13) can be obtained for the case cf 
irrotational flow (-a^ = w) by solving the Laplace equation for 
the velocity potential with an ellipsoidal bounda 3 :^ condition 
(Basset 1961). Also, for -STS O; we get the rigid flow (v = 0), 
Ihus, eq. (13) satisfies all ilie flow requirements. It \'*all be 
shorn -later that the velocity field given by eq, ( 13 ) is unique, 
and JRi corresponds to the angular frequency of the fluid 
particle and has got important significance. 1 


lEhe velocity field defined by eq. (13) is used to obtain; 
the intrinsic energy (H-j^) and the intrinsic angular momen'tum ; 
(1) using eq, (10) with v in place of u. The tr/o are given by i 


^b 


= (?/5)M 2 

i,3 




(a?+a?) 



(Ha) 



( 1 4o) 


- -(4/5) [3 M 


Viihere i, 3 , k are taken in cyclic order. Use of the ahoTe 
equations in eqs. (11) and (12) gives 


H = 

3 

S (Aj^ w? + 

n. a. ^ .X. -u 

(15a) 

= 

vdiere = 

3 

2 (A*, w? + + g! w.-a-. ) 

.>j -LX XX XXX 

■a[ , , 

-r = ^Ri/2; h = -F = - V2 

(15b) 

c. = 

1 

0[/(2Igi) = - 4 . Ijy, a| a^CaJ + 

(15c) 

a;id = 

(IH/ 5 ) (a? + a|). 


Both eqs. ( 15 a) 

and (15b) are in the standard form. 

In both 


the expressions the first term corresponds to pure rotations, 
1he second to the intrinsic particle motion and the third to 
the rotation-intrinsic particle coupling. Clearly, for.A= 0 , 
the flow becomes rigid flow and eqs. (I5a) and. (I5b) reduce to 


= iSi I "ri 4 


t2 _ V t2 „2 
1=1 


and for JfVs w , the flow becomes irrotational aid eqs. 
and (15b) reduce to 


(15d) 

(15e) 

(15a) 



- -(4/5) 


( 1 ‘ 


vAiere i, 3 , i: are talcen in cyclic order. Uvse of the abOTe 
equations in eqs. (11) and (12) gives 


H = 

2 (A- w? + B-ct? + G. w.J^.) 

1 1 1 1 1 1 1 

(15a) 

II 

CM 

2 (a1 w? + BlJiS + c! w.-^. ) 

2.=i ±. X X i XXX 

YiT , ^ , 

(15b) 

Tshere = 

-r = W2! h= -f = -V2 


=1 = 

oycaini) = -+ Igi a| a^(a| + 4)^ 

(15c) 

and = 

(M/5) (a? + a|). 


Both eqs. (15a) 

and (15b) are in the standard form. 

In both 


the expressions the first term corresponds to pure rotations, 
the second to the intrinsic particle motion and the third to 
the rotation-intrinsic particle coupling. Clearly, for-A = 0, 
the flow becomes rigid flow and eqs. (15a) and (I5b) reduce to 


% = i?, I % *1 


(15d) 


2 ^ 

14 = S I 

^ i=1 


2 ^2 
m % 


(I5e) 


and for -Jft= w , the flow becomes irrotational and eqs, (15a) 
and (I5b) reduce to 



% = , 1 ^ I 4 


(15f) 


j2 _ ^ j2 2 


(15g) 


^mere is given by eq* (2) (Bohr and Mottelson 1955) 


^1*2.2 Uniqueness of the Blow 

The uniqueness of the proposed flow given by eq. (15) 
can be demonstrated in a very simple way. 


Let v’ and v” be t\70 solutions satisfying all the 
floT7 requirements, the field W = v’ - v” satisfies the 
conditions ^ • W = ^xW = 0 everywhere, and W*n = 0 on 
the surface, which leads to w = v’ - v” = 0 and hence v is 
unique. 

'^^•2.3 Frequency of Intrinsic Motion 

¥e v/ill now show that the quantity occurring in 
eq. (15) corresponds to the frequency of the fluid element in 
motion. 

B?iting Ca|/(a? + a|) as ^ in cq. (15), we 

obtain v =-»/U x r which means is the angular frequency ^ 

of the fluid element v/ith respect to B and hence is a, 
measure of the frequency of intrinsic motion, 

i 

The relative magnitudes of and w will toll us ' 

whether the system is performing an adiabatic or a non-adiabatic 
motion, lEliis will be discussed later. 



ROTATIOHAL EAMTLTOITIAJr OF. THE SISM 


[Hie deformed even— even nuclei in the rare-earth region 
show a characteristic rotational spectrum, [Therefore, to a 
good approximation, it should be possible to write their 
Hamiltonian in the foim 

5 p 

H = 2 ^i/2% i (16) 

i ""1 • 

vdiere are the components of the effective moment of ineri 

[The Hamiltonian and the angular momentum of the ellipsoi 
with the rotational flow given by eq, (15) are, in their most 
general form, given by eqs, (15a) and (15b) respectively. 

[Hie Hamiltonian as given by eq, (15a) cany in general, 
be put in the rotational foim given by eq, (16) onl.y if it is 
assumed that 


-^i = ^1^1 ('1'^) 

vdiere k^^'s are dimensionless numbers called vorticity para- 
meter, If all k^’s are equal to zero, eq, (13) gives the 
rigid body flow and if they are all equal to unity, v/e get the 
irrotational flow. 

Since Wj^*s are the components of the angular velocity 
of the ellipsoid as a whole, eqs* (13) and (17) imply that the 
intrinsic motion as observed from B is directly coupled with 


the rotation of the ellipsoid as a #iole, which can be considcroc 




to be adiabatic compared to the intrinsic motion only if 
vorticity parameters are large oo335>ared to imity. 

Using the form given by eq. (17) for v/e get the 

rotational Hamiltonian in the form given by eq, (16) with the 
effective moment of inertia given by 

...(18) 

^diere bj^ = a^id is called the shape parameter of the 

ellipsoid. STote that both and eq, (13) go over to the 

known values of the rigid body and the irrotational flow, as 
goes to zero and unity respectively, 

Ihus, the kinetic energy of rotations and the form of 
rotational Hamiltonian compare closely with that of the rigid 
rotator with different expressions for the principal moments 
of inertia for a deformed nucleus. (The Ii^’s occurring in 
eq. (16) are the projections of the total angular momentinn 1 
along the principal axes, and in the operator form they obey 
the commutation relations ^j ^k* Hence, the 

Hamiltonian is quantized in the angular momentum representation 
Ihus, the eigenfunctions of the Hamiltonian operator are nothin 
but the eigenfunctions of the angular momentum operators #iich 
are the usual D-functions , the properties of which are given i 
in Appendix II .4. 



In tile next section we apply eqs. (16) and (18) for the 
^tudy of rotations of a syminetrio and an asynnaetric nuclear 
liquid drop separately. 

TEE CASE OP A SUMBTBIC MD AIT ASYIMEIHIO ITaGIiBAR. .MOP 
Tile aTailable energy and quadrupole moment data for even- 
even nuclei in the rare-earth region indicate that these nuclei 
can "be assumed to "be axially symmetric in their ground state, 
and the effect of asymmetry sets on gradually at higher excita- 
tions. Hence j for the study of these nuclei, it is essential 
to consider the case of a symmetric and an asymmetric drop 
separately* 

The rotational Hamiltonian given "by eq« (16) is exactly 
the same as the rigid rotator Hamiltonian (Chapter III)» except 
that the quantity occurs in place of and hence all 

the results obtained in the case of the rigid rotator must also 
he valid in the present case, 

VI, 4,1 Symmetric Huclear Drop 

If we assume the deformed nucleus to he an axially 
symmetric (about the 3~axis in B) ellipsoidal drop, then 
a^ = ap < a^ and hence Ig-j = Ij |2 > ^^5 ’ since 

+ Ii| = , eq. (16) reduces to 

Hg = (Ii^/2I^) + (E^/2) (1/Ijj3 “ 


(19a) 



’■ ■■’''.cTe K is the pro;5ection of 2j aJong the symmetry axis 
cf the nucleus. She corresponding eigenvalues are given hy 


Eg = I(L+m2/( ) + (k2/2) (I/I1.3 - VIbi) 

ihe values of parameters and occurring in the above 

could be obtained using the knoMm values Of energies of states 
chacacterized by various values of L and K.. It is interest- 
_ng to know that in the case of rotation of a nucleus about an 
axis perpendicular to the syimetry axis (K=0), eqs. <19) will 
contain only the first tessa* Ihis corresponds to the so-called 
ground state rotational band and will be discussed in detail 
in section VI,8,1, 


71,4,2 Asymmetric Ifuclear IXcop ' 

In this case, a^ =|= ag a^. Supposing a^ > a 2 > » 

iien > 1^2 ^ 'possible to obtain the eigen-^ 

values of the rotational Hamiltonian given by eq, (16) exactly. 
However, we obtain them by treating Ihe asymmetry as a small 
perturbation over Idle symmetric state (Davidson 1968), in v/hich 
case eq, (16) can be written in the form 


H = Hg + (1/1^2 - 1/Ig,) (20a) 

Dor ^ 1^2 » "fche last term in eq. (20a) can be treated as 

.. perturbation over the symmetry. Hence, 


< \^ l \ > = - ^2^/2 


<20b) 



Sq, (2013) combined with eq, (20a) will give the eigenvalues 
of H for an slightly asyimietric drop and the eigenfunctions 
are a linear comhination of the eigenfunctions of the symmetri 
drop, 

71.5 ROIAIIOITAI. PLOW IH TWO DIMMSIOIS 

It will he interesting to study the rotational flow in 
two dimensions, as it is much simpler to understand and as it 
gives rise to the rotational spectrum in the ground state 
rotational hand (K=0). 

71,5.1 Iwo-Dimensional Plow as a Particular Case of the 
Plow in Ihree Dimensions 

Consider a tsrpical deformed nucleus of prolate spheroidal 
shape with x^-axis as the axis of symmetry (a^ = a .2 < a^), and ; 
x^-axis as the axis of rotation. OQie flow represented hy eqs. ( 
and (17) becomes two-dimensional flow ^en the fluid elements 
move in planes perpendicular to the axis of rotation ishich 
essentially means ^2 = = 0, Ihis leads to 



with “ ~^2 ^3 / ^^3 ^2^ 

Y&xere a^ and a .2 are respectively the semi -major and the 
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semi-minor axes of tlie prolate spheroid. It is clear from 
eq. (21) that only one vorticity parameter (k^ ) occurs in the 
velocity field. Also, eq. (22) indicates that the flow lines 
are concentric ellipses in planes perpendiciilar to the axis 
of rotation. Oltiis will "be discussed later. 


In the following, we give an alternative method of 
obtaining the velocity field given by eq, (21). 


YI,5*2 Alternative Method 

It 7/ill be shown now tbat starting from anisotropic 
harmonic oscillator forces obtained by including non- rigidity 
in the rigid body forces, it is possible to derive the velocity 
field given by eq. (21). 0 

Supposing the nuclear fluid rotated rigidly about the 
x^-a:cis, the con^onents of the force P acting on a fluid 
element of volume dv taken along the axes in B, are 



f dv w! 


1 ^2 ’ 


P- 


Xr 


f dv w^ 






...(23) 


If, on the other hand, the fluid rotated non-rigidly, 
the velocity v of a fluid element measured in B is not 
necessarily zero. In this case, we assume that the force P 
on dv to be 



dv X2 




ay wf ^3 


p = 0 

..,(24) 
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Sqs. (24) are identical to eqs. (25) for = ^2 ~ Also, 
eqs. (24) give the acceleration of a fluid element in S> if 
S is chosen to coincide with B instantaneously « 

Bor a steady flow of an inconrpressihle fluid, we have 

(av/9t) = 0; f.v = y.u=0* 

Hence, (dV^’fc) = (v*^ )v (25) 

The velocity field v in B satisfies the following equations 
obtained by the usual transfonaation of acceleration in S 
into B (eq, (6) of Chapter V) 

(v.^r )V2 = Xg + 2w^ v^ + w^ ^2 

(26) 

(v.^^ )Vj = ^2 - 2w^ Vg + W^ X5 

and v^ =0. These equations have to be solved to obtain Vg 
and v^ explicitly. 

One of the simplest solutions of eqs. (18) satisfying 
the incompressibility condition is 

Vg = px^ J v^ = qx2 (27) 

ISie Tallies of ’tho coefficienijs p and obtained using 
eq. (27) in eq. (26); and comparing eq, (27) with eq. (22), are 

p = w^[i-Kvi^-'i2)/4 - h n,- ■^ 1 -' V /{] 

= 2a| /(a| + a^) 


(28a) 





i~( n2)/4- H n-j- s( %) a 


2a,^ w^/(a2 3-^) 


(28b) 


alhe algebraic sign before the radica,! is fixed by the rigid 
floT.' reqiiirement, namely, for ^ = "^2 ~ Jai^st have 

A ~ A ~ ^ hence p and g must vanish. Solving for 

and yi^ we get 

Hj = 1 - 4 k^ a.pUl + a^) + 4 k^ a^/Ca^ + 

(29) 

= 1-4 k^a2/(a2 + a|) +4 k^ a^ e|/(a2 + a|)^. 

Using eqs. (28) in eq. (27) and the resultant errprossion in 
eq. (4) we get 


X u = 2|[ 1 - (p - q)/w^ Jw = 2(1 - k^)w, (50) 

(a) Rigid glow ? In this case = ■''I 2 1 '^2 ~ ''^3 ~ 

"^hich leads to p = q = 0 and hence eq. (22) gives the vorticity 
for the rigid flow. 

(b) Irrotational Plow : Using eqs. (28) in oq. (30) we have 

v X u =(^2) ^ ~( n ^ ^ + &C% + ^>2) ‘ 


For the flow to be irrotational the quantity on the right hand 
side must be zero and hence either or must bo negative. 



't'lc shape paranie'fcers of tlie prolaije spheroid o;:te,x 2 ied Iron 
the quadi*upole moment data and iiie valne of 1^, obtained using 
the energy of 1=2 state in eq, (32), 

The rotating non-rigid nucleus alv/ays undergoes c 
centrifugal stretching and hence the corresponding change in 
the potential energy has to be added to H gi'/en by eq[, (o'! a) 
to get the total Hamiltonian of the system, The field exper- 
ienced by a fluid element, represented by the force components 
given by eq, (24), is a conservative field since it satisfies 
= 0. Hence, the force must be derivable from a scalar 
potential field, Eius, the total potential energy’' of the 
system obtained usiiig the force components given by cq. (24), is 

Y = (M wVlO) ( + ^2 a^) (52a) 

vhere K is the mass of the nucleus and and arc 

explicitly given by eqs. (29) in terms of ag, and 

w = It will be shown in 71.8.1 that the potential energy 

of the nucleus as given by eq, (32a) has the familiar Cp /2 
form, 3hus, if the drop is assumed to undergo imiform rota,tions 
with negligible centrifugal stretching, 7 remains a constant 
and hence H given by eq, (31a) can as ?/ell be taken as the 
•total Hamiltonian of the system, 

71.6 SUUPY 01 7jlRIATI0H OP THB MODEL MOMT OP IITERTIA (Ij.) 

She model moment of inertia (1^2.^ given by eq, (13) 
depends on the vorticity parameter k^ and the Gha,pe parameter 



f 


(c) Rotational glow i Por all values k-j , except 0 and 1, 
we obtaip. the rotational flow. But the experimental data 
indicate that 


and hence to obtain rotational flows vhich reproduce the 
experimental values of moments of inertia we must have 0 < < 1. 

VI.5.3 Study of Biiies of Blow 

The flow pattern is studied by plotting the velocity 
con^jonents given by eq. (21) with respect to BCx^XgX^) and 
is shown in Eigure Tl.la. Ohe intrinsic velocity v is zero 
at the origin and increases gradually as we go away from the 
origin and attains its maximum value at the boundary. As the 
direction of velocity field v (clockwise) is opposed to the 
direction of angular velocity w (anti-clockwise) of the whole 
spheroid, each fluid element in motion along concentric elliptic 
paths will experience a certain amount of drag which, in effect, 
results in a sort of slippage of the spheroid, ©lus, v 
effectively decreases w from its rigid body value and hence 
causes a reduction in the rotational kinetic energy of the system 

To study the flow line pattern in S, we have plotted, 
the velocity field u in S(m) and is shown in Rigure Tl.lb. 
The ellipse undergoes a conpression along one of the diameters 
of its director circle followed by a corresponding extension 



Fig.VI.1: Lines of flow in the space (b) 
and body (od fixed coordinate 
systems. 


along the other. Such a motion effectively introduces a diange 
in the orientation of the ellipse in space, [Ehis is true for 
all other ellipses drawn in planes parallel to YZ-plane, ’Ihus, 
the spheroid undergoes a net 2 ?otation in space about the Z-axis. 

YI,5»'4 Rotational Hamiltonian in !I!wo Dimensions 

Ehe constant vorticity rotational flow in two dimensions 
is characterized by = 0 and the vorticity parameter 

for rotations about the x^-axis#- In this case, the rota- 
tional Hamiltonian reduces to the following form if we use ^ 
eqs, (17). (15b) and (I5c) in eq. (15a) 

H = lV(2Ig) (31a) 


vdiere 


-E 


= I 


4- a^a^(a|-a|)^ 

L Q (a^+a2)^-4 k ^ -I 


and the energy eigenvalues are given by 


(31) 


E(I,0) = L(I+l)ilV(2Ij,) 


(32) 


(In the above eq.uations we have dropped the subscript 1 for 
ihe vorticity constant), 

Ihus, the rotational states in this case are characterized 
by K = 0 and L = 0, 2, 4j etc,, giving rise to the ground 
state rotational band, Ihe value of the vorticity parameter (k) 
occurring in eq, (31) can be obtained using the known values of 




It v;ill "be cf theoretical i]atere.-t to 


h . = ( a • /a-. ) « 
variation of with 

for a given k. . 

YI,6,1 Yariation of 


for a given , 


) Vvdth k^ for 0 , given 


and 


h 


udv ti: 


with 


i 


h- 

1 


She variation of ^ function of k^ , for 

various values of h^ 4 1 , is shown in figure VI. 2. lor Id- - 1 
(l^/lj^) = 1 for all values of k^. Also (I^./I^^) gees to 
zero for those values of k. for which 

\ "bi) (33) 


is satisfied. It is clear from the graph that the values of 
k^ corresponding to the minima in the curve goes on increasing 
f 2 ?om 1 as we go from = 1 to h^ = 0. She observed rota- 
tional spectrum of even— even nuclei in the rare-earth region 
indicate that the value of (Ij-iAili^ 0.55. She 

quadrupole moment data indicate that the range of the shape 
parameter h^ is 0.7 to 0.9. In this range of h^ , the values 
of the quantity less than 0.55 are obtained for 

0,8 < k^ < 1,9. However, for ^ ’^he effective moment of 
inertia (Ijji) gees to zero for a finite deformation, which is 
contrary to what is observed. We, therefore, consider the valid 
range of k^ for the nuclei in the rare-earth region as 

0.8 < kj; < 1. 



YI*6.2 Yariation of ) v/ith for a :^iven 

nrne variation of (^ 2 j_Agj_) sis a function of the shape 
parameter h. for various values of k. is sliov.’n in Pig. YI,3. 
I\Fote that for = 1 , the effective moment of inertia hr,s the 
rigid body value for rotational flow, and it lias the value sero 
for the irrotational flow. However, small but finite values of 
effective moment of inertia are obtained for rotational flow 
with kji^ and b^^ slightly less than unity. Such small but 
finite values can also be obtained with the irrotational flow 
but with deformations vdiich are much larger than those for the 
rotational flow* Further, the zero value of the effective 
moment of inertia can also be obtained with rotational flow, 
if varies wdth b^ such that k^^ = 1 when bj^ = 1 , 


■VI. 7 NEED FOR SHE IHTK) HUOJIOH OP BLLIPSOim SHjgES 

Eb.e deformed even-even nuclei in the rare-earth region 
show a rotational spectrum v/ith band structure vdiich closely 
resembles that expected of a spheroid. It is therefore assumed 
that the departure from spheroidal symmetry, if any, is small. 


lo a first order in perturbation, the eigenvalues ol the 
Hamiltonian given by eq. (16) are then given by eq. (20a). 
Since the expected departure from sphericity is small, wc 


consider eq, (20a) in the approximation (1/1^, - 1/I^u) -* 0. 


Ihe values of and 


can then bo obta,inod from the 


energy level data of L=2, E=0(in the ground state band) 



\\ 

\\ 
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kj=:0 


-X X 
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kj = 1.0 Irrotational 
kj = 0.0 Rigid Flow 


V”ki =0.8 

\ * 


[— kj = .99 


\ I « 

vV V-'!/ 


kj = 1.0 


Ilg* Ti;‘3 Taa?iati©a of .(I^^/Ip-) as a ftme-feion of tlie shape 
parameter (hi) ^ for a giTen Talne of the 
Tortieity ooastaat (ki)*. 



and L = K = 2 sdatea usiiid ea* (19^). 3]ne Talnes of 1 -^, 3 , 
obtained are small but finite as shown in Table 71,1. These 
small but finite values of imply that b- = a^j/a^ is 

not unioy as can be seen from Eigure 71.3. Tlie observed omall 
values of can be reproduced by introducing ellipsoidal 


asymmetry, i.e,, b^ ^ 1., The asymmetry is large for = 1 
as compared ti 
Figure 71,2). 


as compared to values of k^ slightly less than unity (see 


Thus, if the higher rotational bands are assumed to be 
precGssional bands, then according to tlie continuum mo col ‘.lie 
nucleus must have ellipsoidal deformation. 


71.8 GAICULATIOITS , 

It is clear from eqs, (19b), (20a) and (20b) that all the 
states in a deformed even-even nucleus can be a-ccounted for 
in terms of either simple rotational motion of a spheroids,! 
drop or in terms of the precessional motion of sji ellipsoidal 
drop. Explicit calculations have been done for tho rotational 
states in the ground state rotation?!, band and in tho oxcited 
state rotational band (the so-called p-vibrationa! bond). Both 
these bands are characterised by K = 0, 

'^1*8,1 Ground State Rotational Band (K = O) 

Rotational states in this bond arc assumed to arise ou.t 
of rotation of the spheroidal nucleus s,bout an axis perpoiidicular 



gABlE 71.1 


Experimental Values of Energies of the L=2 state in the Ground 
State Band (K=0) and L=2 State in the K=2 Band, and the Corres- 
ponding Values of Moments of Inertia 


lucleus 

E(2.0) 

(KeV) 

^E1 

(10“'’^ KeV sec.^) 

E(2.2)* 

(KeV) 

^E3 

(10"''%eV seo.^ 


121.8 

95.9 

1084.0 

7.18 


82.1 

142.4 

1440.0 

5.41 

'5^(}d 

123.1 

94.9 

996.3 

7.82 


89.0 

131.3 

1154.1 

6.75 

15®(}d 

79.5 

146.9 

1187.1 

6,56 

15V 

137.6 

84.9 

890.6 

8,74 

15V 

99.0 

118.0 

946.0 

8.23 

1®>]Dy 

86.8 

134.6 

966.1 

8.06 

16V 

80.7 

144.8 

888.2 

8.77 

15V 

73.4 

159 .'2 

761*8 

10.22 

16V 

102.0 

114.5 

901.0 

8.64 

16V 

91.4 

127.8 

860.3 

9.05 

156ir 

80.6 

145.0 

787.0 

9.89 

168a 

79.8 

146.4 

821.2 

9.48 

ITOet 

79.5 

147.3 

930.0 

8.37 

168yt 

88.0 

132.7 

986.0 

’ -7.90 

i™rb 

84.3 

138.6 

1138.0 

6.84 

1’^V 

78.7 

148.4 

. '-1466^0 

5-31 



TmJE TI > 1 ( Conclude d) 



82,1 

142.3 

1270.0 

6,13 

17% 

88.4 

132.2 

1341,2 

5.81 

178jjj 

93.2 

125.4 

1174.8 

6.63 

CD 

r- 

93.3 

125.2 

1200.2 

6,49 

1®°W 

103.6 

• 112,8 

828,0 

9.41 

1®2W 

100,1 

116.7 

1221.0 

6,38 

1®+W 

111.2 

105.0 

904.0 

8.61 

186 ^ 

122.5 

95.4 

737.0 

10,57 

1®«03 

137.2 

85.1 

767.4 • 

10.15 

1®®0s 

155.0 

75.4 

633.0 

12.13 

19°0s 

180.7 

62.6 

557.9 

13.96 


*Mitsuo Saiai (1970) 



"bo •fclie symmetry’' axis and hence are characterized "by K = 0 
and L = 0, 2, 4» etc. Thus, only one vorticity parameter 
appears in the calculations. The energy eigemralues for 
rotations about the x^-^aris are given by eq, (32) and the 
corresponding moments of inertia are given by eq, (31). 

lor 32 even-even rare-earth nuclei the 2^ state energies 
and the quadrupole moments are knovai (Mariscotti et. al. 1969). 
Assimiing the centrifugal stretching to be small, the former 
data are used to calculate the value of from eq. (32), 

and the latter is used to calculate the values of the parameters 
ag and a^ of the sph'^^’oid using eq. (4) of Chapter II and the 
volume constancy condition. These are then used to calculate 
the value of k from eq. (31). The values of 1^^^ , a2 , a^ , 
P and k for various nuclei in the rare-earth region are given 
in Table VI, 2. Tlie vorticity parameter k is found to vary 
linearly with p from nucleus to nucleus as shown in 
Pigure VI, 4, The empirical relationship foimd from this graph 
is 

k = 1 - 1.5 (34) 

The above expression ■''its to + 10 %• Notice -that for 
spherical nuclei (3 -* O) k is "unity which corresponds to 
the irrotational flow' vvith no rotational spectrum. Eqs, (5) 
of Chapter II v/ith = 0 and eq, (34) are used in eq, (51) 



TABLE VI. 2 


Values of Ie 353 )t * ^2 * ^5 ’ ^ ^ Various Even- 

Even Ifuclei with Spheroidal Shape in the Eare-Earth Region. 


Nucleus 

^expt ^2 

KeV sec^) (fm) 

^3 

(fin) 


k 

152s»l 

95.9 

5.871 

7.620 

0.290 

0.889 

154sm 

142.4 

5.827 

7.857 

0.333 

0,826 

^5+(ja 

94.9 

5.898 

7.651 

0.290 

0.892 


131.3 

5.869 

7.826 

0.523 

0.847 

i5Saa 

146.9 

5.867 

7.953 

0.540 

0,828 


155.1 

5.878 

8*005 

0.348 

0.820 

15Sy 

84.9 

5.934 

7.654 

0*285 

0.907 

158By 

118*0 

5.914 

7.805 

0*310 

0*869 

■'®E5r 

134.6 

5.959 

7i859 

0.310 

0.851 


144.8 

5.952 

7*902 

0.317 

0,842 


159,2 

5.956 

7*990 

0.530 

0.825 


114.5 

5.975 

7.841 

0.303 

0*878 


127.8 

5.989 

7.902 

0,310 

0.865 


145.0 

5.990 

7.994 

0.323 

0,848 

’'®®Er 

146.4 

6.017 

8,018 

0.322 

0.850 

1T0a, 

147,3 

6.057 

8.007 

0.31^V 

0.852 

168^ 

132.7 

6,048 

7.936 

0.303 

0.865 


158.6 

6.065 

7.986 

0*308 

0,862 




TABIiB 71 « 2 (Concluded) 


'•'STb 

148.4 

6.079 

8.044 

0.313 

0,855 

174 ^ 

152.7 ' 

6.119 

8.030 

0-,303 

0.854 

176jt, 

142,3 

6.157 

8.021 

0.295 

0.869 

174 Hf 

128.5 

6,152 

7,943 

0.284 

0*881 

176 jjj 

132.2 

6.173 

7.981 

0.286 

0.879 

178 Hf 

125*4 

6.239 

7.903 

0,262 

0.891 

180 jjj 

125.2 

6.268 

7.916 

0.258 

0.893 

180 ^ 

112.8 

6,206 

7.872 

0.249 

0.904 

182^ 

116.7 

6.324 

7.863 

0.240 

0.905 

1S4, 

105.0 

6.376 

7,822 

0.225 

0.917 

186 ^ 

95.4 

6.407 

7.830 

0.220 

0,926 


85.1 

6.444 

7.740 

0.201 

0,935 

1880s 

75.4 

6.491 

7.710 

0.188 

0.944 

19°0s 

62.6 

6,530 

7,700 

0.180 

0.954 






( 35 ) 


to -eet I™ OTirelr s-- ^ fvjietxon of p ; 

'■“* lii 

= 5Bj^ + 0(p) 1 

•wSiere is tiie mass parameter given "by 

= [] (15 + 8%)^/60Q%'2M Rq . 

!Eh.e potential energy due to centrifugal stretching 
obtained using eqs, (5) of Chapter II with 1^ = 0, and eqs. (29; 

O 

and (34) in eq, (32a) is approximately of the same form (CpV2) 
as that given hy Bohr and Mottelson (1953) to a, lirst oraer on 
p. Thus, the total energy of the system v/ith respect to ohe 
ground state of iiie nucleus is given hy (Eder 1970) 

Ed.O) = IZKl + 3 + - Po)72 (36) 

where is the equilibrium deformation obtained by setting 

( 9B(L}O)/0p ) = 0 for each L. C is the usual stiffness 
parameter obtained by fitting the 4"^ state energy data. These 
values of p^ and C are used to calculate the energies of 
higher rotational states for fifty nuclei in the region 
150 < A < 190. The results are tabulated in Table '7I,3» The 
maximum departure of the calculated value from the ezqjerimental 
value (Mariscotti et. al. 1969) usually occurred for the highest 
IniGwn level. The maximum departure (^E)inax ^ 

for 23 nuclei, less than 2 % for 8 nuclei and more than 2 /q 
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for 19 n-aclei. Ehe values of stiffness constant that we obtain 
are found to be comparable to those of Mosel and Greiner (1968) 
and are considerably greater than those obtained by Mariscotti 
et* al, ( 1969 ). Just as the case of the two mass-point model, 
the stiffness G increases rapidly with the addition of 
neutron pairs and decreases gradually with the addition of prot . 
pairs# 

YI«8,2 Plrst Excited Rotational Band (g -Vibrational Band) 

Rotational states belonging to this band are also 
characterized by K = 0 and L = 0, 2, 4, etc., but they are 
hi^er up in energy with respect to the ground state. Hence, 
in the calculations we have used eg., (36) and the corresponding 
energy minimization condition obtained by putting (8E/8p) * 0 
for each 1. As there is no quadrupole moment data available 
for states belonging to this rotational band, the dependence 
of vorticity parameter k on p was taken to be the same as 
given by eg., (34). 

In the actual calculations we have used the experimental 
values of 1(0"^ ) and E(2'^ ) respectively to obtain the 
stiffness (G) of the harmonic potential, and equilibrium 
def033nation p , for the band by exactly fitting the experi- 
mental energies. Using these values of the parameters, the 
energies of higher rotational states belonging to this band 
have been calculated, !Ehe results, for twenty nuclei in the 
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rare— earijli region for wliicii flie experiiaen'ual values of 3(C ) 

and 3(2"^’) are availalDle, are tabulated in Sable 71. 4. Slie 
rredicted energies for nuclei in tlie middle o.^ xne region are 
found to be in excellent agreement with, the experimental values 
Hot ver37‘ good agreement is obtained for the nuclei in the 
transition region and for neutron deficient nuclei, implying 
thereby that these nuclei do not undergo pure rotations. 

+i 

However, slight changes in the fitting of the ii)(0 ) and 

E(2’^’) values are found to change the predicted energies of 
higher rotational states noticeably, and hence a least squares 
fitting would lead to better results. 


The scarcity of data in this band did not enable us to 
draw any conclusion about the systematics of the stiffness 
constant C as a function of the neutron number end the 
proton number. 


71.8,3 Other Higher Rotational Bands and the Det erminatign 

of the ^ape Parameters of the Huclouis of Ellipsoidal 

Shape 

{[he type of calculations given in the laso two sections 
can be carried on for the rotational states belonging to other 
hig heu rotational bands (K =j= C) assuming them xo arise oux of 
the precissional motion of the nucleus of ellipsoidal shape. 

To do this one has to loiow the values of for values of 

b. = a^/ai^ cyclic). Except for K = 2 band, the 

experimental data for other K =|= 0 bands is very scarce. 



TABL3 VI. 4 


Calculation of the Energies of the Rotational levels in the 
Eirst Excited p -vibrational Band of leformed even-even ITuclei 

in the Rare-Earth Region 



1^ 

0+' 

+* 

2 
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6^’ 

-ut 
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10 

SCO-"’ ) 
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ITucleus 







(Ke?) 

(Mg7 

150m 

1 : 

0.0 

165.0 

465.0 






2 i 

0.0 

165.0 

455.1 

820,3 

1243 

1714 

675.0 

26. 


3 s 

0.2050 

0.2442 

0.2857 

0.3220 

0.3542 

0.3840 



15°siii 

1 s 

0.0 

305.7 

708.9 






2 s 

0.0 

305.7 

725.7 

1208 

1738 

2310 

740.4 

31 .i 


3 s 

0.1200 

0.1900 

0 . 2378 

0.2762 

0.3092 

0.3384 



152sm 

1 ; 

0,0 

126.0 

341 ,0 






2: 

0.0 

126.0 

369. 6 

691.6 

1074 

1507 

685.0 

27. ^ 


3: 

0.2438 

0.2700 

0.3042 

0.3365 

0.3661 

0.3932 



15+sa 

1 : 

0,0 

77,8 

271.1 






2$ 

0.0 
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254.1 

519.3 
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73.: 


3 s 
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Is 

0.0 

79.7 

248.2 





92,11 

2 : 

0.0 

79.7 

260,8 

534.0 

889.7 

1318 

1050.0 


5 s 

0.3353 

0.3391 

0.3472 

0.3579 

0.3702 

0.5832 



15V 

1 1 

0.0 

153.1 

412.3 



1641 

675.6 


2 i 

0.0 

153.1 

428.7 

779.1 

1186 

28,21 

3 s 

0,2089 

0.2435 

0.2823 

0.3168 

0,3475 

0.3754 
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1 s 

0.0 

92.0 

287.0 
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62. 1< 
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gij3LB Vl.4- (Ooncluded) 



1 ! 

0,0 

68.0 

223.9 
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> * 
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174 fj ^ 
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482.0 

805.0 

( 1185 ) 
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3 s 
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tin -aie alDove, 1 and 2 denote tlie experimental and calculated values, 
respectively of the energies E(I+’ )-E(0+* ) and E . in ICeV and 
3 denotes the value of the deformation parameter p , , . Values ’ 
within the parentheses involve a lot of experimental ^ uncertainty 
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Eovrevsr, no .q^Liadrupolo moment data is a-vailaole for these 
states. Hence, v/e adopt a slightly different method to obtain 
the values of as a function of in a self-consistent 

\7e.y, Por this, it is necessary to estimate the valties of 
shape parameters » ^2 ^3* 

In order to obtain b^ , b2 ajid b^ we follow an 
iterative process starting with initial values obtained from 
the limiting case of (l/Igg “ Vlg-]) 0. In this case, the 
values of 1^^ and I^^ can be obtained as mentioned in YI.7 
(refer to Table VI.1), Assuming that the nucler^r volume is 
given by (1,2)\ fm^ and using the quadrupole moment data, the 
shape parameters of the spheroid are obtained. YiTith such values 
of b^ for various nuclei, values of k^ v;ere determined to 
reproduce values of 1^,^ using eq. (31). 'Vaiaos of k^ and 
b^ so obtained could be fitted Tvith the expression 

k^ = 1 - a (1 - b^)2 (37) 

■viSiere i = 1 and a v/as found to be ec^ual to C.71 to vfithin 
10 % • Hote that eq, (37) when expressed in terms of p is 
of the same form as oq. (34). 'Kie values of b^ to fit the 
experimental values of I^,^ are obtained assimiing eq, (37) to 
be valid for i = 3 and using eq. (31). 

!I!hus, by considering the limiting case with ( 1 ^ ’ 

we have obtained the relationship between the vorticity constant 
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and the shape parameter, and an estimate of the value of h^. 

We now consider the case v^here (1/1^2 “ ^ small 

but finite. For an ellipsoid, the quadrupole moment can be 
escpressed in terms of any two shape parameters such as b. , b^ 
or bg , b^ using eq. (6) of Chapter II and the volume 
constancy condition, in the following form 

(Ve = 0.29 Z (A b^/bj)^/^ Q (2 Ai “ 1^3 - 1) + T6 (b^ - 1) ] (3aa) 

Q/e = 0.29 Z (A [ (2b|-l/b^-1 ) + f6 (1 - 1/b^) J (38b) 

Ihus, from the estimate of b^ , estimates of b^ and b 2 are 

obtained using eqs. (38a) and (38b). In the iteration, these 
values of b^ , b 2 » b^ and the parameter a = 0.71 are used 
as the initial values. 

In the case 1^,^ = 1^,2 , the energy E(2,0) of the 
ground state band is a function of I^^ and I |,2 according to 
Eq. (20a). Using eq. (31) and eq. (37), it caai be expressed as 
a function of b^ , b 2 and a. Using estimates of b^ and b 2 
a new value of a is obtained to fit the observed E(2,0) 
value. Ihe energy E(2,2) is a function of 1^,^ , 1^2 and 
Ijj^ according to eq. (20a), Using eqs. (31) and (38) with the 
new value of a , it can be expressed as a function of b^ , b 2 » 
bj. Using the estimates of b^ and b 2 , a ne?/ value of b^ 

is obtained to fit the value of E(2,2). With the new value of 

b^ , new values of b^ and b 2 are found using eqs. (38). 
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We liaYs thus obtained nevr estimates of , h.-, , .uea a. 

Tlie process v;as repeated until self-consistency- v/as achieved. 

Sue procedure was carried o'at for 29 ircclei in the rare-earth 
region for vdiich the enerpgy of K = L = 2 state v;as Inaovm. 
file values obtained are tabulated in Table VI. 5, 

■The self-consistent values of b^ ,bp,b^ and a can 

be tised to calcul?,te the energies of various states in the J:={=G 
bands . 

Thus 5 the IffiM can also be applied for studying rota- 
tional sta,tes in call the encited state rotational bcands. however, 
the data is very scarce for bands with K=}=:0, except for K=2. 

■^1*9 ggSULTS iJ'TD TBEIR GOffl>ARlSOIT WITIi O'THSR IJODHLS 

The energies of rot<aticnal states belonging to the 
ground state and the K=0 excited state rotational bands 
C P —vibra oional band), ca.lculated on the basis of the rotational 
flow model are in guite good agreement with the experiment'.!.! 
values as compared to other models. Since slight changes in 
the fitting of lov/est levels introduces significant changes 
in the energy of higher states, a least squares fitting ^;oulu 
give better results. However, exact fitting procedure gives 
proper systematics of the vorticity?- constant and the stiffness 
parameter. Also, unlike other phenomenological and empirical ' 
models (Diamond, et. al. 1964, Iferiscotti et. al. 1969, 

Sood 1968, Grupta 1969), in the present model, the ImovmL correlatxox 
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TiJBLS YI,5 

Yalues of Shape Parame-fcers , hg , and 
the Parameter a for Ellipsoidal Shape 


Shape Parameters 


mcxeus 

^1 

^2 

^3 

a 


0.759 

1.404 

0.938 

0.584 

15+sm 

0.745 

1 .400 

0.959 

0.829 

15^Gd 

0.759 

1.411 

0.934 

0.560 


0.745 

1.418 

0.948 

0.720 

158Gd 

0.732 

1.439 

0.950 

0.748 


0.764 

1.415 

0.925 

0.498 

158^ 

0.752 

1.415 

0.940 

0.653 


0.746 

1,420 

0.944 

0.713 


0.745 

1.423 

0.942 

0.756 


0.742 

1.428 

0.943 

0.810 


0.754 

1.418 

0.936 

0, 606 


0.745 

1.432 

0.937 

0.636 

'^^Er 

0.743 

1.432 

0.940 

0.711 

168Er 

0.745 

1,426 

0.942 

0.710 


0.748 

1.413 

0.946 

0.715 


0.756 

1,395 

0.948 

0.687 

ITO^b 

0.754 

1.392 

0.952 

0.699 


0.752 

1,388 

0.959 

0.730 



TiBLB YI,3 (Concluded) 
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0.765 

1.366 

0.957 

0.727 

176Hf 

0,769 

1.357 

0.958 

0.690 

178ijf 

0.783 

1.337 

0.955 

0.701 

180gj 

0.784 

1.334 

0.956 

0,691 

180^ 

0,796 

1.328 

0.946 

0.671 


0*800 

1.306 

0.957 

0.710 

184^ 

0.805 

1.311 

0.947 

0.649 

186^ 

0.808 

1.321 

0.937 

0.587 

^860s 

0,820 

1.300 

0.938 

0.5 S 2 


0,826 

1.306 

0.928 

0.531 

'9°0s 

0.818 

1.343 

0.910 

0.412 
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■between the quadru.pole ::ior-ent snd the enengp le”3l data 'las 
been talcen into account in a c^ncistent ’.vay. Ihe waliies of 
stiffness constant needed to fit the rotational states in the 
ground state rotational band e,re cori5)arable to those oltainea 
'bj Mosel and G-reiner (1968) and the Tariation of the same with 
the neutron number and the proton number has approximately if '12 
same trend as those obtained on the basic of other models, 
lihe the two macs-point model, the Tariable moment of inertia 
model, etc. The difference in the Talue of stiffness in the 
ground state band and the K=0 excited state band may be du 2 
to the difference in energies of the L=2 state and the 
corresponding 1=0 state in the respective bands. Further, 
it is ra,ther interesting to note that a flov/ v/ith constant 
vorticity gives a quadratic dependence of the vorticity on 
the deformation parameter ^ or the shape parameter bj^ 
such that for spherical nuclei we obtain irrotatioiial flow 
with no rotational spectrum. The self~conaistGrit calculations 
give the universal dependence of vorticity on the shape para- 
meter although it fits only to -within 20 9o’ 

YI.10 DISGUSSIOh 

The constant vorticity flow given by eq. (13) is unique 
and hence the corresponding Hamiltonian H given by eq. (15a) 
is also unique. In fact, the Hamiltonian in its most general 
form contains a pure rotational energy term, an intrinsic particle 
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energy term emd the cov.plinc tetvreezi the tv^c, it go€Z OTor 
to the rigid body case (Hj^) for* ^ = 0 aiii the irro-fational 
flow case (Hp) for ?/, Ihe relative magnitudes of 

and w determine the strengiii of rotation- intrinsic particle 
coupling. The assumption that the frequency of intrinsic 
particle motion is directly proportional, to the rotational 
frequency (refer to eq, (17)) in order to put the general 
Hamiltonian H given by eq, (15a) into the faiuiliar rotaiionai 
form given by eq, (16), leads to the conclusion that the 
intrinsic motion is directly coupled to the rotatioral motion 
meaning thereby tha.t the intrinsic motion always accompanies 
rotational motion or vice~versa, !I!he vorticity parameters 
(k^) vdiich are introduced as proportionality constants play 
an important role in identifying the type of motion. Thus, 
when k^'s are coc^iarable to unity, the motion is non-adiabatic 
and becomes adiabatic only if kj^' s are ver3r large compared to 
unity. However, for k^ > 1 the effective moment of inertia 
Ijji goes to zero for a finite deformation vAiich is contrary 
to what is observed. Thus, adiabatic notion seems to be 
unphysical in case of a rotational flow with constant vorticity. 


Another important outcome of the RHIvI is tha,t small but 
finite values of effective moment of inertia are obtainable 
with kj: and b^ slightly less than unity (refer to figures 
VI. 3 and VI. 4). Thus, it is possible to account for the observed 
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values of moments of inertia eJong the priiicipal axes of 
the ellipsoid, on the basis of the BMl, Althougii, in case of 
a spheroidal nucleus, the model gives the rigid bo(3y moment 
of inertia along the x^-axis v/hich remains to be explained, 
introduction of slight asymmetry about the x*-axis decreases 
the moment of inertia considerably from its rigid body value* 
On the other hand, the irrotational flow model gives a 2 ero 
value of moment of inertia along the x^-axis of a spheroidal 
nucleus, and requires relatively large asymmetry along the 
x^-axis to get the same value of moment of inertia. Thus, 
we conclude that a spheroidal nucleus can undergo rotations 
about its sjnmnetry axis giving rise to a finite value of 
moment of inertia, and in the case of a flow ^idth constant 
vortlcity, such rotations give rise to the rigid body moment 
of inertia* Ihe observed small but finite values of moment 
of inertia about the sjnmetry axis can be reproduced 

by introducing slight ellipsoidal asymmetry, i*e,, b^ =|= 1. 
However, the zero value of moment of inertia for spherical 
nuclei is always obtained since the variation of with 

b^ is such that ^ = 1 vdien b^^ = 1 , 

Further, it is clear from eq* (35) that for axially 
symmetric nuclei rotating about an axis peipeiidicular to the 
symmetry axis, the REM. moment of inertia depends quadratically 
on the deformation parameter p . Okie same result was obtained 



"by Bohr and Mottelson (1953) lor "the case of irirotational 
flow. Hov/ever, in case of the EKI the terns in higher powers 
of p make a significant contribution to the value of 

XJ 

In fact, a comoarison of the mass parameter (B,-,) obtained on 
the basis of the 2 ?otational flow model \vith that of the 
hydrodynamic mass parameter (Bp) gives (3g/%) =7.1. fhus, 
ti:ie REM has been successful in increasing the value of the 
hydrodynamic mass parameter, thereby implying that the 
effective mass participating in the rotation with a rotationai 
floy; is about seven times the mass participating in the rota- 
tion with an irrotational flow, 

inother interesting consequence of tlie IiEM is that it 
gives the potential energy due to centrifugal stretching in 
the form (p-p^) which is in accordance with that obtained by 
Bohr and Mottelson (1953) and Diamond et. al. (I364)j but is 
contrary to the (p -p^) dependence of Ivlariscotti et. al. (19G9). 
However, the contribution of the potential energ 3 ^ to the total 
energy is relatively small and hence does not make much of a 
difference in the calculations. Lastly, for tlie general three 
dimensional flow the nature of the effective force and the 
corresponding potential energy function are quite complicated. 

Ihus, the REE'I provides an energy Hamiltonian for axially 
symmetric nuclei purely in terms of its deformation or its shape 
parameter, in the form 
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where 


H(P) = 


-r <■ 
dJ 




n o 
o ^ 


= 5 C (15 + Sit) V 60071 [^1 +0(p)]. 
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GHAP PER - iril 
SUMMARY 

Phe quaniruin meclxanical treatment of the simple two 
mass-point model gires the rotational-vibrational energy of 
a deformed ev^'-even rare-earth nucleus in a closed form 
with an I-dependent zero-point energy term which contributes 
considerably to the energy of hi^ spin states with no break- 
down in the rotational structure. Olie model successfully 
correlates the intrinsic quadrupole moment data with the 
moment of inertia more satisfactorily than the rigid spheroid 
or the hydrodynamic models, Phe two mass-point mass parameter 
(Bj) is about half of the hydrodynamic value and the 

corresponding moment of inertia (l^) is about a third of the 
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■<!!> 

rigid spheroid -value tox values oi However, 

■unlike does not go to zero with p, Further, the 

model gives proper values of intrinsic quadrupole moments 
and B(E2; 4"*" "*• 2'^)/B(E2| 2^ -* o"^) for transitions of the 
type Ij + 2 -* L* (Ehe variation of the stiffness (0) of the 
harmonic potential with the neutron number and with the 
proton number are in accordance with those found on the basis 
of various other models. 

Use of the enharmonic Morse potential v;ith exactly 
the same stiffness as that of the' harmonic potential did not 
give any noticeable improvement in the fitting of energy levels. 

Both the potentials give almost the same values of the 
parameters (1^ and C), and the largest value of the stiffness 
is obtained for ^Hf , ail almost * rigid rotor’ . Both the 
potentials yield excellent results vdien a least squares fitting 
procedure is adopted., 

Ihe states belonging to the first excited p~vibrational 
band can also be analysed on the basis of the tvro mass— point 
model, Howeven, the data available is scarce. 

Ihe t-wo mass-point model fails to account for the 
S-t 3 ^e systematics observed for the nuclei of interest, [Ehe 
main drawback of the model is that it does not- give zero 
moment of inertia for the case of spherical nuclei* However,- 
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this has heeii taken care of in the generalised tvro mass-point 
model with a rotationally in-variant core of variable radius, 
and a quantum mechanical treatment of the same ‘TOuld, perhaps, 
yield improved results as compared to the simple tv/o mass— 
point model, 

!Ehe rotational flow model (EJM) based on the principles 
of continuum mechanics with a constant vorticity flow in three 
dimensions gives satisfactory values of energies and deforma- 
tions for the even— even nuclei in the rare-earth region, !IIhe 
proposed rotational flow in three dimensions is unique and 
hence the corresponding Hamiltonian is also -unique. Ihe 
Hamiltonian (H) and the angular moment-urn (1) thiis obtained, 
in their most general form, contain three terms corresponding 
to pure rotational motion, the intrinsic-particle motion 
and the coupling between the two, and they go over to the 
rigid body case and the irrotational flow case under appropriate 
limits of the vorticity, She general Hamiltonian can be put 
into the familiar rotational form (1 /2I) only if it is 
assumed that the frequency of intrinsic-particle motion is 
directly coupled to the rotational frequency through the 
vorticity parameters. Thus, the velocity field directly 
depends on the rotational frequency of the ellipsoidal bo-undary 
and hence the adiabatic motion seems to be not possible in 


this case. 
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By far the most interesting aspect of the RIM is the 
study of variation of a function of for a 

given hj^ , and as a function of h^ for a given k^. Whereas 
the former gives the physically valid range of vorticity 
constant (k) for the nuclei of interest as 0,8 < k < 1, the 
latter shows that small hut finite values of moment of inertia 
along the ssrmmetry axis could he obtained for values of k^ 
and hj^ slightly less than unity. Both the graphs indicate 
if the higher rotational hands are assumed to he precessional 
hands, then according to the continuum model the nucleus 
must have ellipsoidal deformation. Thus, the rotational states 
in the higher rotational hands can he analysed on the basis of 
the RIM, Ihrther, the graphs indicate that vdthiii the frame- 
work of the RIM, a spheroidal nucleus can undergo rotations 
about its symmetry axis giving rise to a finite vo-lue of 
moment of inertia, and in the case of a flow with constant 
vorticity, such rotations give rise to the rigid body moment 
of inertia. The observed small hut finite values of moment of 
inertia (Ig^) about the symmetry axis can he reproduced by 
introducing slight ellipsoidal asymmetry. However, the zero 
value of moment 'of inertia for spherical nuclei is alwa-ys 
obtained since the variation of k^^ with hj^ is such tha-t 
kj_ = 1 when h^ = 1 . 

Also , for axially symmetric nuclei rotating about an 
axis perpendicular to the symmetry axis, the RIM moment of 
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inertia depends quadrat ically on the deformation parameter p, 
HoweTer-, the terms in higher powers of p iiiahe a significant 
contrihution to the value of I. In fact, the RHl mass 
parameter, thus obtained, is about seven times the hydrodynamic 
mass parameter implying thereby that the effective mass 
participating in the rotational motion with rotational flow 
is about seven times the mass participating in the rotational 
motion with an irrotational flov/. 

In the case of a simple two dimensional flow, the BIM. 
is found to give a p^ dependent potential energy term arising 
'Out of the centrifugal stretching of the rotating nucleus, 

Ihus, the RBI provides an energy Hamiltonian for axially 
s3rmmetric nuclei purely in terms of its deformation or its 
shape parameter, 

Ihe RM fails to account for the S-type systematics 
observed for the nuclei of interest and an attempt to reproduce 
the same is desirable. Application of the Rill to the study of 
deformed nuclei in other regions of the periodic table may be 
of interest while examining the validity of the model in general. 
As an extension of the work presented in this thesis, it will 
be interesting to study the nature of the effective force and 
the corresponding potential energy function in the case of 
the general three dimensional constant vorticity flow. IPurilier, 

as the deformed nuclei can not only rotate but also vibrate, it 

<1 



will lie of use to investigate the possibility of including 
shape vibrations of the ellipsoidal drop by considering its 
shape parameters to be time dependent. Shis vhll be useful 
to account for the spectra of spherical nuclei on the basis 
of the rotational flov;- model. 
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APPENDIX I 


It is ImoTOi that the energy of rotational states in the 
gronnd state rotational hand increases with the angular momentum 
I of the state in accordance with the 1(1+1 ) rule. In this 
case the moment of inertia is defined as 



dE 

d tL(L+1)J 



( 1 . 1 ) 


and the corresponding rotational frequency is defined as 


■B. w = 


dE 


( 1 . 2 ) 


dVflT'ra)"' 

In order to investigate the variation of (2 I/Rl^) with ^^w^, 
the quantities on the ri^t hand side of eqs. (I.l) and (1,2). 
are obtained by replacing the derivative vfith the coi*re spending 
differential quotient evaluated between the spin values D and 
Ii-2, Thus, using 


dE = E(L) - E(D~2) 

d(D(L+l)) = D(D+1) - (I--2) (D-O = 2(2L-1) 


we obtain 

^ r e(D) - E(I>-2) KeV^ 

(2D-1)^ L -I 

and 

21 ^ 2(2D-1) ■ • 

E(L) - E(L-.2) 

vfeich are eqs. (8) and (9) of Chapter II. 
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iPPMPIX - II 

^•1 Po obtain the azigular momenirum components in their operator 
form 

We have the correspondence relation 

( ) < — > ( <{) 0 f ) <— “> ( ) 

So one can write 

a^- = ^ 0 f ) 

Taking the time derivative 



Wow, hy comparing these w . * s with those given hy eq. (3a), 

O 


we get 

9a^ 


Sin tJf ; 

o>| o> 

CD 

I) 

Cos i|r 1 

9a^ 

w “ ° 

Sin 9 

9oc Q 


Cos llf 5 

9a2 

-Sin "f 1 

9a o 

W" ~ 

Sin 9 

w- = 

11 

o 

da- 

5? = 

Cos 9 

9aj 

= 0 ; 

9a ■2 

D 4 

• 


Writing <}) = ^ a 2 ) yields 

9 9a-j 9 9a2 9 9oc^ 9 

9^ 9<{) 9a^ 9^ 9a2 9^ 9a^ 

Using the above equations, we get 
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9 — n4,^/-N 9 


^ = Sin e Sin T}! + Sin 9 Cos f + Cos 9 ^ 


Similarly one can obtain 


^ = Cos f Sin * ^ 




5^ 


and 


3 

Bf 


3 

Bo^ 


Solve for (3/3a.) to get eqs, (19a)* Similarly, we can write 
J 


a 


j (<j) 0 f) 


Taking the. time derivative 


3a 


“X = 


= IT 't' 1e^ ® + 5|- ^ 


3 a 


Y * 

r\ 


3a, 


ETow, comparing these j ^ with those given by eq.s. (31*) 
we have 


3a-CT* 

ST 

3 tty 

w 

Sag 

w 


= 0 


= 0 


3ax 

BeT 

3ay 

B9“ 

Da^ 


Cos ({) 
Sin (j) 


^“x 

FfT 

3ay 

Ff 

9a rr 


= 0 


= Sin 9 Sin ({) 
= -Sin 9 Cos (|) 
= Cos 9 « 


Eollowing the same procedure as before we get eqs. (I9b)* 


Consider 


■2 t2 . t2 . t2 _ -2 

j = + 1j2 + ^-T = 


^82 a2 82 




Using eqs, (I9a) we get 

1 2 rtj.^2i4r 4 2 ,1 a 2 


3*^ _ Sin'^f 3“^ , 2* 3^^ , „^,2 ^ 3 2Sin f Cos f Cot 

— — = -v i t-x + Cos f — sr + Sin f Cot 9 -yw - qT 7 ; '~ 5 1 

aaf Sin^s ar 3^ ^ 


2 

+ Sin^f Got 9 5 ^ + Sin f Cos f L l±Q°| ,. 9 . 1 ^ 


Sin^ 9 


Similarly, 


4 = 2°§!1 ^ + aosSf cot2e -4 + gg - ^ - J -. g?? .. | - g - ° - ^ - § 


aag sin^e aij)^ 


w 


30^ 3f 

‘ + Cos^Tjr Got-0 ^ - Sin f Cos f L l- tS2|4§I ^ 

^2 


Sin 0 


Sin^0 


Bf 


and 


3a j 


32 

iF 


Adding the three leads to 

( 20 ), 


Sin 0 3f' 


Similarly, working with a^, ay* we will end up with the 


same result for 1“^. 


II* 2 Commutation Relations 

Using the operator form of and ly, Ly, Lg 

the following relations can easily he verified 


and 

Ifz’hl = 


and cyclic in 1,2,3 


and cyclic in Z,Y,Z, 


(21) 


The above relations are used to verify 


n2 , 1 


. t" -r -r “Hh* 
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* T^fcere the subscript •+’ indicates that the quantity is an 
anti-commutator, 

Prom the ahoTe it follov/s that 

= 0 

and 

= 0 

Also the above relations can be used to verify 

= ° . 

[l2,l|] = p.llj = ^ _ 


(21a) 


and 


(211:)) 


1^*3 n^o Obtain the Rotation Operator , R((t) 6 T|J) Explicitly 
In order to obtain R(^ 0 explicitly, it is 
necessars’- to study the effect of rotation on a given function 
f(XIZ), Consider a rotation of a point fixed in BCx-jX^x^) 
frame through an angle a about the common Z-x^ axes. Then 

X-| = X Cos a + Y Sin a 

Xg = “X Sin a + X Cos a 


x^ — Z , 





^ X - 


Row R(a) f(XX2}) = fCx^XgX^) = f(X Cos a+Y Sin a,-Z Sin a+Y Cos 

Por infinitesimal rotation a is small and hence a Taylor 
expan*sion of f can be made about a = 0, irehich is much easier 
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if vie work in spherical polar coordinates (r 0 4 ) , vdaere 

X = r Sin 9 Cos <j) 
y = r Sin 9 Sin f , 

Z = r Cos 0 


-a 


Then R f(xyz) = R f(r, 0,4)= f(r,e, 4 ~a) = e f(re 4 ) 

•vdiere the exponential is i^ie well-known symbolic form of the 
infinite Tayloj* series. 

If Ii is the orbital angular momentum then 


= 


15. 3 

r ^ 


Ihus 


R(a) f(XrZ) = 


T , 

e ^ ^ f(xrz) 


If e^ is the unit vector along the Z-axis, then 


R(a) f 


f 


vdiere ijr is the wave function describing the system, !Ehis can 
be written in a more general form 


R(a) f = 


~^(5*n) 


where R(a) ssrmbolizes a rotation through an angle a about 
any axis and n is a unit vector along that axis. 

If a is small, vie can expand the above equation 
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= [ 1 - ^ (1-5) + 0(<x2) ] f 

giving tto first order change in "Ij due to R(a) as 
- ^ (L,n)ljj. Pbr this reason the three angular momentum operators 
Lj, 1 y» also referred to as the infinitesimal rotation 

operators » (L.n) will simply he some linear combination of 

these operators, the coefficients depending on the direction 
of n of the axis of rotation* Thus, it is possible to write 


the individual rotations 

J^(<I^)>r(9)j R(t) explicitly. 



R((i)) = 

„ n. Z 
e 

R(9) = 

e * 


-ill 

r(t;o = 

e ^ 3 


so that 

R((t) 0 f) = R(f) rL(0) R(<()) 






A simple geometrical argument (Rose, 


R((t) 0 f) 



1957) leads to 


i© 



1 


Z 


which is equation (25). 


Il,4 Pro r qrties of 3>-functions 

Here, we will list some important properties of the 
D-functic.s without going into the details of derivation. 
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(a) Taking the complex conjugate of eq.. (25) of the text 




= <1^1 e^ ^ ^ ^ 1IE> 


= (f(l!vl) 

^ i(h _ 
= ( 






^(ei^ + fi2) 


^ T|r(LE)J (Change of notatio; 


Tlr(lM), e 


iM»(|) 


= e 


(IMj e 


I 

■=5r J 


ILK)' 


iE]! 


t(LK) y^( Since 
J is uni 






unitary) 




iM(j) 


m' 


d:;;.(e) 


iKf 


-ie j. 

Tflhere d^^(e) = (LMi e ^ ^ jLK) 


(25a) 

(25h) 


By group theoretical methods Wigaer has shown that 

3^(6) = f '('M)";' ii.+r)r'u-i')'!Tir-r)T 

(_1)P (COS |)a+E-M-2p (s^ |,M-Ef2p 

, J : 

p pl(l-M-p)I (L+K-p)r Tp+M-K)! 

T/diere the summation extends over all integral values of p for 
which the arguments of the factorials are non-negative, Notice 
that d^(6) is real, that is why the complex conjugation in 
eq_. (2513) has been omitted. Thus, our definition of <3 ^ (8-! 
is in acc'^rdance with definition given hy Rose, although our 
1 ^( 3 ? ) is the complex conjugate of I ^ ( V ) given hy Rose 
(1957). 
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(b) Properties of d^(0) 

Since the inverse of rotation expC-iGli^j-) is expCiGIi^-), 
a rotation about X throng -G means 

T/\4iich leads to 

'i^(®) = (-1)^^ *^-1,1, -K^®^ (25c) 

Por the particular case of 0 = tc we have the very simple result 
djg-(it) = (-1)^“^ 6(M+K, 0) (25d) 

In terms of I ^ ( oJ ) we have 

34(t))* = (-1)“-^ (250) 

The rotation matrices possess the usual orthonormality p’^^'^oerties 

(25f) 

Coupling rules for the are given by 

= I (l^L2M^M2i™)(:^L2i;^K2llK) ( 25g) 

vdiere It = 1^+L2 » ^ ~ M^+M 2 and Z = E^+K 2 * 


Using orthogonality relations we can invert (25g) to give 
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= 2 


MiKi 


(!■ 




. . . ( 25I1) 


In eq» (24) of the text, the functions ijl will in genera, '’onsists 
of a rad-'-al and angular part which are simply the sphericc,! harmonics 
Yjjg, In ''he case of pure rotations tie radial part does not change 
and hence can he dropped out of the equation (24). Ihus, we get 




(25J) 


where ( 0 } (j) ) indicate the spherical polar coordinates the 
S-frame o"d ( 0%9') in the B-frame, 

Integrals of products of B-functions over an unit sphere are 
given hy 




2 ^. 


^ Ig) S(K.|K2) 
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and 
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= o(m^,M2+m^) 6(k-,,e:2+k^) ,'^'1^.^ 

...(25k) 

Bor eith'- ’ M or K = 0 the ) reduce to spherical 

harmonic 
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11,5 lo determine tlie action the relabeling transformations 

HI. 5 *5 have upon these functions, we use 


(a) Since the operation is the rotation % about the 

2- axis 

= < 1^3 = 0 1^2 “ ^ 
and using '’^’0) “ (-1)^"^^ %,-E 

we get 

) ~ ^ (25ii) 

p 

(b) Since T| is two successive rotations of n/2 about the 

3- axis clearly 


D^(0 0 %) 


-±Z% 
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and hence 

o T * iK75 T * 

^2 ^ = e ^ (25p) 



